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Abstract 



Recent equations of state for dense nuclear matter are discussed with possible 
phase transitions arising in neutron stars such as pion, kaon and hyperon konden- 
sation, superfluidity and quark matter. Specifically, we treat the nuclear to quark 
£f~) ' matter phase transition, the possible mixed phase and its structure. A number of 

numerical calculations of rotating neutron stars with and without phase transitions 



are given and compared to observed masses, radii, temperatures and glitches. 
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1 INTRODUCTION 



1.1 The past, present and future of neutron stars 



The discovery of the neutron by Chadwick in 1932 prompted Landau to predict 
the existence of neutron stars. The birth of such stars in supernovae explosions 
was suggested by Baade and Zwicky 1934. First theoretical neutron star cal- 
culations were performed by Tolman, Oppenheimer and Volkoff in 1939 and 
Wheeler in 1960. Bell and Hewish were the first to discover a neutron star 
in 1967 as a radio pulsar. The discovery of the rapidly rotating Crab pulsar 
in the remnant of the Crab supernova observed by the Chinese in 1054 A.D. 
confirmed the link to supernovae. Radio pulsars are rapidly rotating with pe- 
riods in the range 0.033 s < P < 4.0 s. They are believed to be powered by 
rotational energy loss and are rapidly spinning down with period derivatives 
of order P ~ 10~ 12 — 10~ 16 s. Their high magnetic field B leads to dipole 
magnetic braking radiation proportional to the magnetic field squared. One 
estimates magnetic fields of the order of B ~ 10 11 — 10 13 G. The total number 



of pulsars discovered sofar has just exceeded 1000 ahead of the turn of the 
millenium and the number is increasing rapidly. 

A distinct subclass of radio pulsars are millisecond pulsars with periods be- 
tween 1.56 ms< P < 100 ms. The period derivatives are very small corre- 
sponding to very small magnetic fields B ~ 10 8 — 10 10 G. They are believed to 
be recycled pulsars, i.e. old pulsars with low magnetic fields that have been 
spun up by accretion preserving their low magnetic field and therefore only 
slowly spinning down. About 20 - almost half of the millisecond pulsars - are 
found in binaries where the companion is either a white dwarf or a neutron 
star. Six double neutron stars are known sofar including the Hulse- Taylor PSR 
1913+16. The first binary pulsar was found by Hulse and Taylor in 1973 and 
by measuring the inward spiralling one could determine all parameters in the 
binary system as both masses, orbital periods and period derivatives, orbital 
distances and inclination. Hereby general relativity could be tested to an un- 
precedented accuracy. The binary neutron stars all have masses in the narrow 
interval 1.3 — 1.5M , which may either be due to the creation process or that 
heavier neutron stars are unstable. 

With X-ray detectors on board satellites since 1971 almost two hundred X- 
ray pulsars and bursters have been found of which the orbital period has been 
determined for about sixty. The X-ray pulsars and bursters are believed to 
be accreting neutrons stars from high (M^1OM ) and low mass (M^,1.2M ) 
companions respectively. The X-ray pulses are most probably due to strong ac- 
cretion on the magnetic poles emitting X-ray (as northern lights) with orbital 
frequency. The X-ray burst are due to slow accretion spreading all over the 
neutron star surface before igniting in a thermonuclear flash. The resulting (ir- 
regular) bursts have periods depending on accretion rates rather than orbital 
periods. Recently, the "missing link" between bursters and pulsars has been 
discovered. It is the low mass X-ray burster XTE J1808-369 where also 401Hz 
X-ray pulsations have been detected [1]. The radiation from X-ray bursters 
is not blackbody and therefore only upper limits on temperatures can be ex- 
tracted from observed luminosities in most cases. Masses are less accurately 
measured than for binary pulsars. We mention recent mass determinations for 
the X-ray pulsar Vela X-l: M = (1.9 ± O.1)M , and the burster Cygnus X-2: 
M = 1.8±O.2)M , which will be discussed below. A subclass of six anomalous 
X-ray pulsars are slowly rotating but rapidly spinning down indicating that 
they are young with enormous magnetic fields, B ~ 10 13 G, and thus named 
"Magnetars" [2]. Recently, quasi-periodic oscillations (QPO) have been found 
in 11 low mass X-ray binaries. The QPO's set strict limits on masses and radii 
of neutrons stars, but if the periodic oscillations arrive from the innermost 
stable orbit [3], it implies definite neutron star masses up to M ~ 2.3M . 

Non-rotating and non-accreting neutron stars are virtually undetectable. With 
the Hubble space telescope one single thermally radiating neutron star has 



been found [4]. Its distance is only 160 pc from Earth and its surface temper- 
ature is T ~ 60 eV. From its luminosity one deduces a radius of the neutron 
star R < 14 km. In our galaxy astrophysicists expect a large abundance ~ 10 8 
of neutron stars. At least as many supernova explosions have occurred since 
Big Bang which are responsible for all heavier elements present in the Universe 
today. The scarcity of neutron stars in the solar neighborhood may be due to 
a high initial velocity (asymmetric "kick") during their birth in supernovae. 
Recently, many neutron stars have been found far away from their super- 
nova remnants. Such "invisible" neutron stars have probably been detected by 
gravitational microlensing experiments. Future microlensing observation will 
determine the population of such dark matter objects in the galactic halo. 

From the view of physicists (and mass extinctionists) supernova explosions are 
unfortunately rare in our and neighboring galaxies. The predicted rate is 1-3 
per century in our galaxy but the most recent one was 1987A in LMC. With 
luck we may observe one in the near future which produces a rapidly rotat- 
ing pulsar. Light curves and neutrino counts will test supernova and neutron 
star models. The rapid spin down may be exploited to test the structure and 
possible phase transitions in the cores of neutron stars [5-7]. 

The recent discovery of afterglow in Gamma Ray Bursters (GRB) allows de- 
termination of the very high redshifts (z > 1) and thus the enormous distance 
and energy output E ~ 10 53 ergs in GRB. Most previous GRB models cannot 
produce these energies. Only neutron star mergers or black holes may be able 
to produce such violent events. In comparison supernovae produce ~ 10 51 erg's 
mainly in neutrinos whereas a long distance missile has (only) ~ 10 31 ergs of 
explosive energy. 

The marvelous discoveries made in the past few decades will continue as nu- 
merous satellite experiments are running at present and will be launched. 
History tells us that the future will bring great surprises and discoveries in 
this field. 



1.2 Physics of neutron stars 



The physics of compact objects like neutron stars offers an intriguing inter- 
play between nuclear processes and astrophysical observables. Neutron stars 
exhibit conditions far from those encountered on earth; typically, expected 
densities p of a neutron star interior are of the order of 10 3 or more times the 
density pd ~ 4- 10 11 g/cm 3 at 'neutron drip', the density at which nuclei begin 
to dissolve and merge together. Thus, the determination of an equation of 
state (EoS) for dense matter is essential to calculations of neutron star prop- 
erties. The EoS determines properties such as the mass range, the mass-radius 
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Fig. 1. Possible structure of a neutron star. 

relationship, the crust thickness and the cooling rate. The same EoS is also 
crucial in calculating the energy released in a supernova explosion. Clearly, 
the relevant degrees of freedom will not be the same in the crust region of a 
neutron star, where the density is much smaller than the saturation density 
of nuclear matter, and in the center of the star, where density is so high that 
models based solely on interacting nucleons are questionable. These features 
are pictorially displayed in Fig. 1. Neutron star models including various so- 
called realistic equations of state result in the following general picture of the 
interior of a neutron star. The surface region, with typical densities p < 10 6 
g/cm 3 , is a region in which temperatures and magnetic fields may affect the 
equation of state. The outer crust for 10 6 g/cm 3 < p < 4 ■ 10 11 g/cm 3 is a 
solid region where a Coulomb lattice of heavy nuclei coexist in /^-equilibrium 
with a relativistic degenerate electron gas. The inner crust for 4 • 10 11 g/cm 3 
< p < 2 ■ 10 14 g/cm 3 consists of a lattice of neutron-rich nuclei together with 
a superfluid neutron gas and an electron gas. The neutron liquid for 2 • 10 14 
g/cm 3 < p < -10 15 g/cm 3 contains mainly superfluid neutrons with a smaller 
concentration of superconducting protons and normal electrons [8] . At higher 
densities, typically 2 — 3 times nuclear matter saturation density, interesting 
phase transitions from a phase with just nucleonic degrees of freedom to quark 
matter may take place [9]. Furthermore, one may have a mixed phase of quark 
and nuclear matter [5,10], kaon [11] or pion condensates [12,13], hyperonic 
matter [5,14-19], strong magnetic fields in young stars [20,21] etc. 

The first aim of this work is therefore to attempt at a review of various ap- 
proaches to the equation of state for dense neutron star matter relevant for 
stars which have achieved thermal equilibrium. Various approaches to the EoS 
and phases which may occur in a neutron star are discussed in section 2 while 
an overview of the thermodynamical properties of the mixed phase and pos- 
sible phases in neutron stars are presented in sections 3 and 4. 



Our second aim is to discuss the relation between the EoS and various neu- 



tron star observables when a phase transition in the interior of the star occurs. 
Astronomical observations leading to global neutron star parameters such as 
the total mass, radius, or moment of inertia, are important since they are sen- 
sitive to microscopic model calculations. The mass, together with the moment 
of inertia, are also the gross structural parameters of a neutron star which 
are most accesible to observation. It is the mass which controls the gravita- 
tional interaction of the star with other systems such as a binary companion. 
The moment of inertia controls the energy stored in rotation and thereby the 
energy available to the pulsar emission mechanism. Determining the possible 
ranges of neutron star is not only important in constraining the EoS, but has 
important theoretical consequences for the observational prediction of black 
holes in the universe. Examples are the galactic black hole candidates Cyg 
X-l [22] and LMC X-3 [23], which are massive X-ray binaries. Their mass 
functions (0.25 M and 2.3 M ) are, however, smaller than for some low-mass 
X-ray binaries like A0620-00 [24] and V404 Cyg [25] , which make better black 
hole candidates with mass functions in excess of three solar masses. There is 
a maximum mass a non-rotating neutron star can have. There is however no 
upper limit on the mass a black hole can have. If, therefore, one can find a 
dense, highly compact object and can argue that its rotation is slow, and can 
deduce that its mass is greater than the allowed maximum mass allowed to 
non-rotating (or slowly rotating) neutron stars, then one has a candidate for 
a black hole. 

Since neutron stars are objects of highly compressed matter, this means that 
the geometry of space-time is changed considerably form that of a flat space. 
Stellar models must therefore be based on Einstein's theory of general relativ- 
ity. Based on several of the theoretical equations of state and possible phases 
of matter discussed in sections 2, 3 and 4, various properties of non-rotating 
and rotating neutron stars are presented in section 5. The relevant equations 
needed for the study of the structure of a neutron star are summarized in this 
section as well, for both non-rotating and rotating stellar structures. There we 
also discuss the obervational implications when phase transitions occur in the 
interior of the star. In addition to studies of the mass-radius relationship and 
the moment of inertia, we extract also analytical properties of quantities like 
the braking index and the rate of slow down near the critical angular velocity 
where the pressure inside the star just exceeds that needed to make a phase 
transition. The observational properties for first- and second-order phase tran- 
sitions are also discussed. Other properties like glicthes and cooling of stars 
are also discussed in section 5. Summary and perspectives are given in section 
6. 

Finally, we mention several excellent and recent review articles covering var- 
ious aspects of neutron stars properties in the literature addressing the in- 
teresting physics of the neutron star crust [26] , the nuclear equation of state 
[13], hot neutron star matter [27] in connection with protoneutron stars, and 



cooling calculations [28]. However, as previously mentioned, our aim will be 
to focus on the connection between the various possible phases of dense neu- 
tron star matter in chemical equilibrium and the implications of first- and 
second-order phase transitions for various observables. 



2 PHASES OF DENSE MATTER 



Several theoretical approaches to the EoS for the interior of a neutron star 
have been considered. Over the past two decades many authors [9] have con- 
sidered the existence of quark matter in neutron stars. Assuming a first order 
phase transition one has, depending on the equation of states, found either 
complete strange quark matter stars or neutron stars with a core of quark mat- 
ter surrounded by a mantle of nuclear matter and a crust on top. Recently, 
the possibility of a mixed phase of quark and nuclear matter was considered 
[5] and found to be energetically favorable. Including surface and Coulomb en- 
ergies this mixed phase was still found to be favored for reasonable bulk and 
interface properties [10]. The structure of the mixed phase of quark matter 
embedded in nuclear matter with a uniform background of electrons was stud- 
ied and resembles that in the neutron drip region in the crust. Starting from 
the outside, the crust consists of the outer layer, which is a dense solid of neu- 
tron rich nuclei, and the inner layer in which neutrons have dripped and form 
a neutron gas coexisting with the nuclei. The structure of the latter mixed 
phase has recently been calculated in detail [26] and is found to exhibit rod-, 
plate- and bubble-like structures. At nuclear saturation density po — 2.8- 10 14 
g/cm 3 there is only one phase of uniform nuclear matter consisting of mainly 
neutrons, a small fraction of protons and the same amount of electrons to 
achieve charge neutrality. A mixed phase of quark matter (QM) and nuclear 
matter (NM) appears already around a few times nuclear saturation density 
- lower than the phase transition in hybrid stars. In the beginning only few 
droplets of quark matter appear but at higher densities their number increase 
and they merge into: QM rods, QM plates, NM rods, NM bubbles, and finally 
pure QM at very high densities if the neutron stars have not become unstable 
towards gravitational collapse. 

In this section we review various attempts at describing the above possible 
phases of dense neutron star matter. For the part of the neutron star that 
can be described in terms of nucleonic degrees of freedom only, i.e. /3-stable 
matter with protons, neutrons and electrons (and also muons), we will try 
to shed light on recent advances within the frawework of various many-body 
approaches. This review is presented in subsection 2.2. For the more exotic 
states of matter such as hyperonic degrees of freedom we will point to recent 
studies of hyperonic matter in terms of more microscopic models in subsec- 
tion 2.4. The problem however with e.g. hyperonic degrees of freedom is that 



knowledge of the hyperon-nucleon or hyperon-hyperon interactions has not yet 
reached the level of sophistication encountered in the nucleon-nucleon sector. 
Mean fields methods have however been much favoured in studies of hyperonic 
matter. A discussion of pion and kaon condensation will also be presented in 
the two subsequent subsections. Superfluidity is addressed in subsection 2.7. 

In general, we will avoid a discussion of non-relativistic and relativistic mean 
fields methods of relevance for neutron matter studies, mainly since such as- 
pects have been covered in depth in the literature see e.g. Refs. [16,17,29,30]. 
Moreover, as pointed out by Akmal et al. [13], albeit exhibiting valuable tu- 
torial features, the main problem with relativistic mean field methods is that 
they rely on the approximation /ir < 1, with /i the inverse Compton wave 
length of the meson and r the interparticle spacing. For nuclear and neu- 
tron matter densities ranging from saturation density to five times saturation 
density, fir is in the range 1.4 to 0.8 for the pion and 7.8 to 4.7 for vector 
mesons. Clearly, these values are far from being small. The relativistic mean 
field approximation can however be based on effective values for the coupling 
constants, taking thereby into account correlation effects. These coupling con- 
stants have however a density dependence and a more microscopic theory is 
needed to calculate them. 

Our knowledge of quark matter is however limited, and we will resort to phe- 
nomenological models in subsection 2.8 in our description of this phase of mat- 
ter. Typical models are the so-called Bag model [31] or the Color-Dielectric 
model [32]. 

However, before proceeding with the above more specific aspects of neutron 
star matter, we need to introduce some general properties and features which 
will enter our description of dense matter. These are introduced in the first 
subsection. The reader should also note that we will omit a discussion of the 
properties of matter in the crust of the star since this is covered in depth 
by the review of Ravenhall and Pethick [26]. Morever, for neutron stars with 
masses ~ 1.4M or greater, the mass fraction contained in the crust of the 
star is less than about 2%. We will therefore in our final EoS employ results 
from earlier works [33,34] for matter at densities < 0.05 fm -3 . 



2.1 Prerequisites and definitions 



At densities of 0.1 fm" 3 and greater, we will in this work require properties of 
charge neutral uniform matter to be made of mainly neutrons, protons, elec- 
trons and muons in beta equilibrium, although the presence of other baryons 



will be discussed as wellQ. 

In this section we will merely focus on distinct phases of matter, such as 
pure baryonic matter or quark matter. The composition of matter is then 
determined by the requirements of chemical and electrical equilibrium. Fur- 
thermore, we will also consider matter at temperatures much lower than the 
typical Fermi energies. The equilibrium conditions are governed by the weak 
processes (normally referred to as the processes for /^-equilibrium) 

b 1 ^b 2 + l + vi b 2 + 1 -> h + ui, (1) 

where b\ and b 2 refer to e.g. the baryons being a neutron and a proton, re- 
spectively, I is either an electron or a muon and V\ and v\ their respective 
anti-neutrinos and neutrinos. Muons typically appear at a density close to 
nuclear matter saturation density, the latter being 

n w0.16 ±0.02 far 3 , 

with a corresponding binding energy £ for symmetric nuclear matter (SNM) 
at saturation density of 

S = B/A = -15.6 ± 0.2 MeV. 

In this work the energy per baryon £ will always be in units of MeV, while 
the energy density e will be in units of MeVfm -3 and the number density^ n 
in units of fm -3 . The pressure P is defined through the relation 

P = n 2 — = n— -e, (2) 

on on 

with dimension MeVfm -3 . Similarly, the chemical potential for particle species 
i is given by 



& = [-ZT ' ( 3 ) 




with dimension MeV. In our calculations of properties of neutron star matter in 
/3-equilibrium, we will need to calculate the energy per baryon £ for e.g. several 
proton fractions x p , which corresponds to the ratio of protons as compared to 
the total nucleon number (Z/A), defined as 

x P = ^, (4) 

n 



1 In this work we will also set G = c = h = 1, where G is the gravitational 
constant. 

2 We will often loosely just use density in our discussions. 
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where n = n p + n n , the total baryonic density if neutrons and protons are 
the only baryons present. In that case, the total Fermi momentum kp and 
the Fermi momenta kp p , kp n for protons and neutrons are related to the total 
nucleon density n by 

2 

= x p n + (1 — x p )n 
I i 

= 3^ ,+ 3^ v (5) 

The energy per baryon will thus be labelled as S(n, x p ). S(n, 0) will then refer 
to the energy per baryon for pure neutron matter (PNM) while £ (n, |) is the 
corresponding value for SNM. Furthermore, in this work, subscripts n,p,e,fi 
will always refer to neutrons, protons, electrons and muons, respectively. 

Since the mean free path of a neutrino in a neutron star is bigger than the 
typical radius of such a star (~ 10 km), we will throughout assume that 
neutrinos escape freely from the neutron star, see e.g. the work of Prakash et 
al. in Ref. [27] for a discussion on trapped neutrinos. Eq. (1) yields then the 
following conditions for matter in (3 equilibrium with e.g. nucleonic degrees 
freedom only 

Hn = /J>p + He, (6) 

and 

n p = n e , (7) 

where //» and rii refer to the chemical potential and number density in fm~ 3 of 
particle species i. If muons are present as well, we need to modify the equation 
for charge conservation, Eq. (7), to read 

Tip fl e ~r '^, 

and require that /i e = /i M . With more particles present, the equations read 

E (< + O = E fe + "j:) > (8) 

i i 

and 

Hn = hVi + QiHi, (9) 

where bi is the baryon number, q^ the lepton charge and the superscripts (±) 
on number densities n represent particles with positive or negative charge. To 
give an example, it is possible to have baryonic matter with hyperons like A 
and X~' 0,+ and isobars A~ ,0 ' +,++ as well in addition to the nucleonic degrees of 
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freedom. In this case the chemical equilibrium condition of Eq. (9 ) becomes, 
excluding muons, 

fJ'-E- = HA~ = Hn + He, 

HA = HZ = /! A = /i ra , 

HT.+ — HA+ — H~p = Hn — He, 

HA++ = Hn - 2/i e . (10) 



A transition from hadronic to quark matter is expected at high densities. 
The high-density quark matter phase in the interior of neutron stars is also 
described by requiring the system to be locally neutral 

(2/3K " (l/3)n d - (l/3)n s - n e = 0, (11) 

where n u ^,s,e are the densities of the u, d and s quarks and of the electrons 
(eventually muons as well), respectively. Morover, the system must be in (3- 
equilibrium, i.e. the chemical potentials have to satisfy the following equations: 

Hd = Hu + He, (12) 

and 

Hs = Hu + He- (13) 

Eqs. (11)-(13) have to be solved self-consistently together with e.g. the field 
equations for quarks at a fixed density n = n u + n d + n s . In this section we will 
mainly deal with distinct phases of matter, the additional constraints coming 
from the existence of a mixed phase of hadrons and quarks and the related 
thermodynamics will be discussed in Section 3. 

An important ingredient in the discussion of the EoS and the criteria for 
matter in /^-equilibrium is the so-called symmetry energy S(n), defined as the 
difference in energy for symmetric nuclear matter and pure neutron matter 

S{n) =£{n,x p = 0)-£{n,x p = 1/2). (14) 

If we expand the energy per baryon in the case of nucleonic degrees of freedom 
only in the proton concentration x p about the value of the energy for SNM 
(x p = |), we obtain, 

1 1 d 2 £ 
£(n, x p ) = £(n, x p = -) + - -j^ (n) (1 - 2x p ) 2 + ..., (15) 

where the term d 2 £jdx 2 is to be associated with the symmetry energy S(n) in 
the empirical mass formula. If we assume that higher order derivatives in the 
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above expansion are small (we will see examples of this in the next subsection), 
then through the conditions for /5-equilbrium of Eqs. (6) and (7) and Eq. (3) 
we can define the proton fraction by the symmetry energy as 



2 \V3 



he [37r 2 nx p j = 4S(n) (1 - 2x p ) , (16) 

where the electron chemical potential is given by \x e = hckp, i.e. ultrarelativis- 
tic electrons are assumed. Thus, the symmetry energy is of paramount impor- 
tance for studies of neutron star matter in /^-equilibrium. One can extract 
information about the value of the symmetry energy at saturation density n 
from systematic studies of the masses of atomic nuclei. However, these results 
are limited to densities around no and for proton fractions close to |. Typical 
values for S(n) at n are in the range 27 — 38 MeV. For densities greater than 
no it is more difficult to get a reliable information on the symmetry energy, 
and thereby the related proton fraction. We will shed more light on this topic 
in the next subsection. 

Finally, another property of interest in the discussion of the various equations 
of state is the incompressibility modulus K at non-zero pressure 

dP 

K = 9 ^ < 17) 

The sound speed v s depends as well on the density of the nuclear medium 
through the relation 



2 dP dPdn ( K 



cj de dn de \9(m n c 2 + S + P/n) J 

It is important to keep track of the dependence on density of v s since a super- 
luminal behavior can occur at higher densities for most non-relativistic EoS. 
Superluminal behavior would not occur with a fully relativistic theory, and it 
is necessary to gauge the magnitude of the effect it introduces at the higher 
densities. This will be discussed at the end of this section. The adiabatic con- 
stant T can also be extracted from the EoS by 



2.2 Nucleonic degrees of freedom 



A major part of the densities inside neutron stars can be well represented by 
nucleonic degrees of freedom only, namely the inner part of the crust to the 
outer part of the core, i.e. densities ranging from 0.5 to 2 — 3 times nuclear 
matter saturation density, There is a wealth of experimental and theoretical 
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data, see e.g. Ref. [35] for an overview, which lend support to the assumption 
that nucleons do not loose their individuality in dense matter, i.e. that proper- 
ties of the nucleon at such densities are rather close to those of free nucleons. 
The above density range would correspond to internucleon distances of the 
order of ~ 1 fm. At such interparticle distances there is little overlap between 
the various nucleons and we may therefore assume that they still behave as 
individual nucleons and that one can absorb the effects of overlap into the 
two nucleon interaction. The latter, when embedded in a nuclear medium, is 
also different from the free nucleon-nucleon interaction. In the medium there 
are interaction mechanisms which obviously are absent in vacuum. As an ex- 
ample, the one-pion exchange potential is modified in nuclear matter due to 
"softening" of pion degrees of freedom in matter. 

In order to illustrate how the nucleon-nucleon interaction is renormalized in 
a nuclear medium, we will start with the simplest possible many-body ap- 
proach, namely the so-called Brueckner-Hartree-Fock (BHF) approach. This 
is done since the Lippmann-Schwinger equation used to construct the scat- 
tering matrix T, which in turn relates to the phase shifts, is rather similar 
to the G-matrix which enters the BHF approach. The difference resides in 
the introduction of a Pauli-blocking operator in order to prevent scattering to 
intermediate particle states prohibited by the Pauli principle. In addition, the 
single-particle energies of the interacting particle are no longer given by kinetic 
energies only. However, several of the features seen at the level of the scattering 
matrix, pertain to the G-matrix as well. Therefore, if one employs different 
nucleon-nucleon interactions in the calculation of the energy per baryon in 
pure neutron matter with the BHF G-matrix, eventual differences can be re- 
traced at the level of the T-matrix. We will illustrate these aspects in the next 
subsection. More complicated many-body terms and relativistic effects will be 
discussed in subsections 2.2.2 and 2.2.3 



2.2.1 From the NN interaction to the nuclear G-matrix 

The NN interactions we will employ here are the recent models of the Nijmegen 
group [36], the Argonne Vis potential [37] and the charge- dependent Bonn 
interaction (CD-Bonn [38]). In 1993, the Nijmegen group presented a phase- 
shift analysis of all proton-proton and neutron-proton data below 350 MeV 
with a x 2 P er datum of 0.99 for 4301 data entries. The above potentials have 
all been constructed based on these data. The CD-Bonn interaction has a x 2 
per datum of 1.03 and the same is true for the Nijm-I, Nijm-II and Reid93 
potential versions of the Nijmegen group [36]. The new Argonne potential Vis 
[37] has a x 2 P er datum of 1.09. 

Although all these potentials predict almost identical phase shifts, their math- 
ematical structure is quite different. The Argonne potential, the Nijm-II and 
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the Reid93 potentials are non-relativistic potential models defined in terms of 
local potential functions, which are attached to various (non-relativistic) oper- 
ators of the spin, isospin and/or angular momentum operators of the interact- 
ing pair of nucleons. Such approaches to the NN interaction have traditionally 
been quite popular since they are numerically easy to use in configuration 
space calculations. The Nijm-I model is similar to the Nijm-II model, but it 
includes also a p 2 term, see Eq. (13) of Ref. [36], which may be interpreted as 
a non-local contribution to the central force. The CD-Bonn potential is based 
on the relativistic meson-exchange model of Ref. [39] which is non-local and 
cannot be described correctly in terms of local potential functions. 

For a given NN interaction V, the i?-matrix (or ii'-matrix) for free-space two- 
nucleon scattering is obtained from the Lippmann-Schwinger equation, which 
reads in the center-of-mass (cm.) system and in a partial- wave decomposition 

Rf*{kk'u) = V^(kk') + W ^V^{kq)—^—R^f{qk'ul (20) 

with IV and kk! the orbital angular momentum and the linear momentum of 
the relative motion, respectively. T z is the total isospin projection. The angular 
momentum J and total spin S are represented by the variable a. The term 
Hq represents the kinetic energy of the intermediate states. The phase-shifts 
for a given partial wave can be calculated from the on-shell matrix element 

a 2 

of R, which is obtained by setting q = q' = qo with oj = — Q -, m n being the 
mass of the nucleon. Since all of the above interactions reproduce the same 
phase-shifts, the corresponding on-shell matrix elements of R calculated from 
these various potentials are identical as well. However, due to the way the 
potentials are constructed, their off-shell properties may be different. This 
was discussed in detail in Refs. [40,41]. In those works the authors showed 
that especially for the 1 S'o and 3 Si- 3 /}i channels, the CD-Bonn and Nijm-I 
interactions which include the effects of non-localities, yield a more attractive 
interaction for the free scattering case. For D-waves and higher partial waves 
the various potentials were almost equal while there were still differences for 
P- waves. 

We now turn the attention to the application of such NN interactions in a 
nuclear medium. First we will therefore employ a as simple as possible many- 
body scheme, in order to preserve a link between the preceeding discussion on 
the NN interaction and the solution of the Lippmann-Schwinger equation. As 
stated above, it will thus suffice to employ the BHF method. 

Following the conventional many-body approach, we divide the full Hamil- 
tonian H = T + V, with T being the kinetic energy and V the bare NN 
interaction, into an unperturbed part H = T + U and an interacting part 
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Hi = V — U, such that 

H = T + V = H + if 7 , 

where we have introduced an auxiliary single-particle (sp) potential U. If £/ is 
chosen such that Hi becomes small, then perturbative many-body techniques 
can presumably be applied. A serious obstacle to any perturbative treatment 
is the fact that the bare NN interaction V is very large at short inter-nucleonic 
distances, which renders a perturbative approach highly prohibitive. To over- 
come this problem, we introduce the reaction matrix G given by the solution 
of the Bethe-Goldstone equation (in operator form) 

G{u) = V + VQ ^tt^QG, (21) 

where uo is the unperturbed energy of the interacting nucleons and Q is the 
the Pauli operator which prevents scattering into occupied states. The Pauli 
operator is given by 

its k Tm k ~> k Tn 
0, otherwise, 

in the laboratory system, where kp defines the Fermi momenta of the proton 

(tj = 1/2) and neutron (tj = —1/2). For notational economy, we set |k m | = k m . 

The above expression for the Pauli operator is in the laboratory frame. In the 
calculations of the G-matrix, we will employ a Pauli operator in the center- 
of-mass and relative coordinate system. Further, this Pauli operator will be 
given by the so-called angle- average approximation, for details see Ref. [42]. 
Eq. (21) reads then (in a partial wave representation) 

Gf'(kk'Kco) = V,f-(kk') (23) 



I" 



The variable K is the momentum of the center-of-mass motion. Since we are 
going to use an angular average for the Pauli operator, the G-matrix is di- 
agonal in total angular momentum J . Further, the G-matrix is diagonal in 
the center-of-mass orbital momentum L and the total spin S, all three vari- 
ables represented by the index a. The variable a differs therefore from the 
definiton of the i?-matrix, where K — 0. Three different G-matrices have to 
be evaluated, depending on the individual isospins (tit 2 ) of the interacting 
nucleons (||, — \ — \ and —\\)- These quantities are represented by the total 
isospin projection T z in Eq. (23). The different G-matrices originate from the 
discrimination between protons and neutrons in Eq. (22). The term ifo in the 
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denominator of Eq. (23) is the unperturbed energy of the intermediate states 
and depends on k, K and the individual isospin of the interacting particles. 
Only ladder diagrams with intermediate two-particle states are included in 
Eq. (21). The structure of the G-matrix equation in Eq. (23) can then be di- 
rectly compared to the i?-matrix for free NN scattering, Eq. (20). Therefore, 
as discussed below, eventual differences between various potentials in a finite 
medium should be easily retraced to the structure of the i?-matrix. It is also 
obvious that one expects the matrix elements of G to be rather close to those 
of R with only small deviations. These deviations originate from two effects 
which reduce the contributions of second and higher order in V to the G- 
matrix as compared to their contributions to R. One is the above-mentioned 
Pauli quenching effect: the Pauli operator Q in (23) restricts the intermedi- 
ate particle states to states above the Fermi energy. The second one is the 
dispersive effect: the energy denominators in (23) are defined in terms of the 
single-particle energies of nucleons in the medium while the corresponding de- 
nominators of (20) are differences between the energies of free nucleons. Since 
the absolute values for the energy differences between nucleons, which feel 
the mean field of the nuclear system, are larger than the energy differences 
between the kinetic energies, also this dispersive correction reduces the attrac- 
tive contributions of the non-Born terms. As a result, the matrix elements of 
G tend to be less attractive than the corresponding matrix elements of R, see 
e.g. Refs. [40,41] for further details. 

We use a continuous single particle (sp) spectrum advocated by Mahaux et 
al. [43]. It is defined by the self-consistent solution of the following equations: 

k 2 
£ . = ti + Ui = -*- + u h (24) 

2m 

where m is the bare nucleon mass, and 

Ui = Yl (ih\G{E = ei + e h )\ih) AS . (25) 

h<kp 

In Eqs. (24)- (25), the subscripts % and h represent the quantum numbers of 
the single-particle states, such as isospin projections r« and Th, momenta ki 
and kh, etc. The sp kinetic energy is given by ti and similarly the sp potential 
by u^ 

Finally, the non-relativistic energy per nucleon S is formally given as 

1 k 2 

g = 1 V ^ + (26) 

-^ £ (hh'\G(E = e h + e h ,)\hh>) As . 

h<kp,h'<kp 
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In this equation we have suppressed the isospin indices for the Fermi momenta. 
Eq. (26) is actually calculated for various proton fractions x p , and is thereby a 
function of both density n and x p . We will therefore in the following discussion 
always label the energy per particle as £(n,x p ). 

In the limit of pure neutron matter only those partial waves contribute where 
the pair of interacting nucleons is coupled to isospin T = 1. Due to the an- 
tisymmetry of the matrix elements this implies that only partial waves with 
even values for the sum I + S, like l S , 3 P etc. need to be considered in this 
case. For proton fractions different from zero, in particular the case of sym- 
metric nuclear matter, also the other partial waves, like 3 Si — 3 D\ and X P\ 
contribute. In a BHF calculation the kinetic energy is independent of the NN 
interaction chosen. We will then restrict the following discussion to the poten- 
tial energy per nucleon U, the second term in the RHS of Eq. (26). Putting the 
contributions from various channels together^, one obtains the total potential 
energy per nucleon U for symmetric nuclear matter and neutron matter. These 
results are displayed in Figs. 2 and 3 as functions of density n for the the CD- 
Bonn interaction [38] (solid line), the three Nijmegen potentials, Nijm-I (long 
dashes) , Nijm-II (short dashes) and Reid93 (dotted line) [36] and the Argonne 
Vi 8 [37] (dot-dashed line). 

The differences between the various potentials are larger for the energy in nu- 
clear matter. This is mainly due to the importance of the 3 Si- 3 Di contribution 
which is absent in pure neutron matter. This is in line with previous inves- 
tigations, which showed that the predicted binding energy of nuclear matter 
is correlated with the strength of the tensor force, expressed in terms of the 
.D-state probability obtained for the deuteron (see e.g. [42]). This importance 
of the strength of the tensor force is also seen in the calculation of the binding 
energy of the triton in Ref. [38,44]. The CD-Bonn interaction yields a binding 
energy of 8.00 MeV, the Nijm-I potential gives 7.72 MeV while the Nijm-II 
yields 7.62 MeV, the same as does the new Argonne potential [37]. Typically, 
potentials with a smaller D-state probability have a weaker tensor force and 
exhibit therefore a smaller quenching of the non-born terms in Eq. (23). More- 
over, the fact that the CD-Bonn interaction and the Nijm I potential include 
effects of non-localities, yields also a further attraction from the central force 
both in the 3 Si- 3 Di and the singlet 1 Sq channels. The latter explains the addi- 
tional difference in nuclear matter between the Nijm I and the Nijm II, Reid93 
and Argonne Vis potentials as well as part of the difference seens in Fig. 3 for 
pure neutron matter (PNM). For both SNM and PNM there are also addi- 
tional differences arising from P waves, notably for the Argonne potential in 
PNM, where the difference between the Argonne Vis interaction the Reid93 
and the Nijm II potentials is mainly due to more repulsive contributions from 
P-waves. For higher partial waves, the differences are rather small, typically 
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Fig. 2. Potential energy per particle IA for symmetric nuclear matter as function of 
total baryonic density n. 
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Fig. 3. Potential energy per particle IA for pure neutron matter. 

of the order of few per cent. The differences seen in Figs. 2 and 3 should also 
be reflected in the symmetry energy defined in Eq. (14). This is seen in Fig. 4 
where we display the symmetry energy for the above potentials as function of 
density n. From the differences in symmetry energies one would then expect 
that properties like proton fractions in /3-stable matter will be influenced. This 
in turn has important consequences for the composition of matter in a neu- 
tron star and thereby eventual phases present in dense matter. In Fig. 5 we 
display the corresponding proton fractions obtained by calculating the energy 
per particle through the G-matrix of Eq. (23) and imposing the equilibrium 
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Fig. 5. Proton fraction x p for /3-stable matter as function of density n 

conditions of Eq. (6) and including muons, see again Ref. [45] for further de- 
tails. From Fig. 5 one notices that the potential with the largest symmetry 
energy, the CD-Bonn interaction, is also the one which gives the largest proton 
fractions. This means in turn that the so-called direct Urea process can occcur 
at lower densities. For the CD-Bonn interaction this happens at 0.88 fm~ 3 , for 
the Nijm I it starts at 1.25 fm -3 while for the Reid93 interaction one reaches 
the critical density at 1.36 fm -3 . The Argonne potential allows for the direct 
Urea process at a density of 1.05 fm~ 3 . For the Nijm II we were not able to 
get the direct Urea process for densities below 1.5 fm -3 . 
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Fig. 6. Fermi momenta for neutrons, protons, electrons and muons in /3-stable matter 
with the CD-Bonn interaction. 

It is also interesting to notice that the symmetry energy increases rather 
monotonously for all potentials. This means that the higher order deriva- 
tives in Eq. (15) can be neglected and that we can, to a good approximation, 
associate the second derivative d 2 £ /dxi with the symmetry energy S(n) in 
the empirical mass formula. With this mind, one can calculate the proton 
fraction employing the theoretically derived symmetry energy shown in Fig. 4 
using the simple formula of Eq. (16). A good agreement is in general obtained 
with the above simple formula, see e.g. Ref. [45] and the parametrization in 
subsection 2.3. 

With the results for the proton fractions of Fig. 5 in mind, we plot in Fig. 
6 the Fermi momenta for electrons, muons and nucleons obtained with the 
CD-Bonn interaction, the other potentials yield qualitatively similar results 
although the proton fraction is slightly smaller. From this figure one notices 
that muons appear at a density close to the saturation density of nuclear 
matter, as expected. This can easily be seen if one were to calculate /3-stable 
matter with non-interacting particles only, see e.g. Ref. [46]. 

Of further interest is the difference in energy per particle £(n,x p ) for pure 
neutron matter (PNM) and matter in /^-equilibrium. In Fig. 7 we display 
S(n, Xp) for PNM and /3-stable matter for results with the CD-Bonn interaction 
only since the other potentials yield qualitatively similar results. Obviously, as 
seen from Fig. 7 the energy per particle for /5-stable matter yields a softer EoS, 
since the repulsive energy per particle in PNM receives attractive contributions 
from the T z = channel. This is reflected in the corresponding pressure as 
well (lower panel of the same figure) and will in turn result in neutron stars 
with smaller total masses compared with the PNM case. 
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Fig. 7. Upper panel: BHF energy per nucleon for pure neutron matter and /3-stable 
matter obtained with the CD-Bonn interaction. Lower panel : the corresponding 
pressure P. 

We end this subsection by plotting in Fig. 8 the energy density per baryon 
S (including the contribution from leptons) for /3-stable matter. Since the 
proton fractions are not too large, see Fig. 5, the most important contribution 
to e and £ stems from the T z = — 1 channel and the contribution from the 
nuclear tensor force, especially via the 3 Si and 3 Si- 3 Di contributions with 
T z = 0, plays a less significant role than what seen for the potential energy 
in symmetric nuclear matter in Fig. 2. Thus, the main contribution to the 
differences between the various potentials arises from the T z — 1 channel. For 
the Tz = 1 channel, the new potentials yield also results for /3-stable matter 
in close agreement. 
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Fig. 8. The total energy in /3-stable matter as function of density n. 

It is thus gratifying that the new NN interactions yield similar energies per 
particle in neutron star matter. As we will see in the next subsection also, 
more sophisticated many-body calculations at the two-body level yield rather 
similar results for neutron matter, since the strong nuclear tensor force in the 
T z = channel is not present in pure neutron matter or in a less important 
way in /3-stable matter. However, contributions coming from real three-body 
interaction or relativistic effects, may alter this picture. This is the topic of 
the next two subsections. 
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2.2.2 Higher-order many-body calculations 

The previous subsection served to establish the connection between the nucleon- 
nucleon interaction and the simplest many-body approach possible, namely 
the summation of so-called ladder diagrams by the Brueckner-Hartree-Fock 
(BHF) method. We will label these calculations as lowest-order Brueckner 
(LOB) theory. This allowed us to see how the free NN interaction gets modi- 
fied in a nuclear medium. Eventual differences in e.g. the EoS for neutron star 
matter could then be retraced to properties of the various NN interactions. 
The new high-quality NN interactions have also narrowed the differences at 
higher densities observed in the literature for older interaction models, see e.g. 
the discussion in Ref. [40]. However, it is well-known that the BHF method 
in its simplest form is not fully appropriate for a description of dense matter. 
More complicated many-body terms arising from core-polarization effects, ef- 
fective three-body and many-body diagrams and eventually the inclusion of 
three-body forces are expected to be important at densities above no- More- 
over, LOB theory suffers from other pathologies like the lack of conservation 
of number of particles [47]. The need to include e.g. three-body interactions is 
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Fig. 9. First order contribution to the ground-state energy shift AEq in BHF theory. 

seen already at the level of the triton since all the above potentials underbind 
the triton, albeit the discrepancies which existed earlier have been reduced. 
To give an example, the old Reid potential [48] gave a binding energy of —7.35 
MeV while a precursor to the CD-Bonn interaction, the Bonn A interaction 
[39] resulted in —8.35 MeV. The modern potentials yield results in a range 
from -7.6 to -8.0 MeV. 

Here we will therefore review three possible improvements to LOB theory 
and discuss the resulting equations of state in detail. The first improvement 
will be to consider the summation to infinite order of the chain of particle- 
particle hole-hole diagrams (PPHH) [49-52]. Thereafter, we will discuss the 
recent three-hole line results of Baldo and co-workers [53,54] and finally the 
calculations with three-body forces as well by Akmal et al. [12,13]. 

In brief, the summation of PPHH diagrams means that the Pauli operator in 
Eq. (22) is extended in order to prevent scattering into intermediate hole-hole 
states as well. This means that in addition to summing up to infinite order 
diagrams with particle-particle intermediate states, we will now sum, still to 
infinite order, a larger class of diagrams containing hole-hole intermediate 
states as well. 

For LOB, as discussed in the previous subsection, the effective two-body inter- 
action in nuclear matter is given by the G-matrix, which includes all ladder- 
type diagrams with particle-particle intermediate states to infinite order. The 
ground-state energy shift AE in terms of the G-matrix is represented by the 
first-order diagram of Fig. 9 and reads 

AE^ OB = J2n a n b (ab\G(uj = e a + e b )\ab). (27) 

ab 

In Eq. (27) the n's are the unperturbed Fermi-Dirac distribution functions, 
namely nk =1 if k < kp and =0 if k > kF where kp is the Fermi momentum. 

The Brueckner-Hartree-Fock G-matrix contains repeated interactions between 
a pair of "particle" lines, as illustrated by the diagrams of Fig. 10. Note that 
they are so-called Goldstone diagrams, with an explicit time ordering. The 
third-order diagram (a) of Fig. 10 is given by 

,-^. / \ /-*-\3 * abmn* rnnrs*rsab /r>o\ 

Diag.(a) = (-) -, -. (28) 

* ( e a + Q> — t m — e n J(e a + e 6 — e r — e s ) 
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(a) 

Fig. 10. Goldstone diagrams contained in the Brueckner-Hartree-Fock G-matrix. 



(P) 

Fig. 11. Goldstone diagram with repeated interactions between hole lines. 

Here m, n, r, s are all particle lines, and V^i represents the anti-symmetrized 
matrix elements of the NN interaction V. 

To the same order, there is also a diagram with hole-hole interactions, as shown 
by diagram (f3) in Fig. 11, where a, b, r, s are all hole lines. This diagram is 
not included in standard LOB calculations of nuclear matter, for the following 
reason. In earlier times, nuclear-matter calculations were based on, by and 
large, the so-called hole-line-expansion. The essence of the hole-line approach 
is that diagrams with (n+1) hole lines are generally much smaller than those 
with n hole lines. With this criterion, diagram (/?) which has 3 (independent) 
hole lines would be negligible compared with diagram (a) which has 2 hole 
lines. Thus the former could be neglected. To investigate the validity of this 
criterion, it may be useful to actually calculate diagrams like (f3). 

A motivation behind the PPHH-diagram method of nuclear matter is thus to 
include diagrams with hole-hole correlations like diagram (/3) to infinite order 
For further details on how to obtain an effective interaction and the energy 
per particle in neutron star matter, see e.g. Refs. [50,52,55]. 
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Fig. 12. The contributions from the bubble (BUB), ring (RING), and higher-order 
three hole-line diagrams (HIGH) employing the continuous choice. Results are for 
symmetric nuclear matter, employing the Argonne V14 model for the NN interaction. 



The next set of diagrams which can be included is the summation of so-called 
three-hole line diagrams through the solution of the Bethe-Fadeev equations, 
originally pionereed by Day [56] and recently taken up again by Baldo and 
co-workers [53,54]. The whole set of three hole-line diagrams can be grouped 
into three main sets of diagrams. The so-called ring diagrams, the bubble di- 
agram, and the so-called higher-order diagrams, with an arbitrary number of 
particle lines. Each of these three-hole line contributions are quite large, see 
Refs. [53,54] and Fig. 12, but there is a strong degree of cancellation among 
the various terms. Thus, the total three-hole line contribution turns out to be 
substantially smaller than the two-body (i.e. two hole-line) contribution. This 
can be seen in Fig. 12 where we plot the various three hole-line contributions 
and their total contribution, all with the continuous choice for the single par- 
ticle energies. The interesting feature of the calculations of Baldo et al. [54], is 
that the results with the continuous single-particle choice lead to three hole- 
line results which are rather close to the total two hole-line results. Diagrams 
of the PPHH type are however not included. In Ref. [57] these were estimated 
to be of the same sign in both the standard and the continuum choice. In 
spite of this methodological progress in perturbative approaches, there are 
still classes of diagrams which need to be summed up. As discussed by Jack- 
son in Ref. [58] , one can prove that there is a minimal set of diagrams which 
need to be summed up in order get the physics of a many-body system right. 
This set of diagrams is the so-called Parquet class of diagrams [59] where ring 
diagrams and ladder diagrams are summed up to all orders in a self-consistent 
way. Since it is rather hard to sum this set of diagrams in a practical way, 
other many-body methods like the coupled-cluster ansatz [60] and optimized 
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hypernetted-chain theory [61] provide a systematic approximation method to 
sets of Feynman diagrams that cannot be calculated exactly. In the remainder 
of this section will therefore focus on a variational method based on the hy- 
pernetted chain summation techniques developed by Pandharipande, Wiringa 
and co-workers. The approach was developed in the 1970s [62], particularly to 
include the effects of many-body correlations, presumably important in dense 
neutron star matter. Calculations performed since then have confirmed that 
many-body clusters make significant contributions to the binding energies of 
equilibrium nuclear matter and light nuclei [56,62-65]. We will base our discus- 
sion of the EoS on the recent results of Akmal, Pandharipande and Ravenhall 
[12,13,66]. We will refer the reader to the latter references and Ref. [62] for 
more details. In brief, the variational wavefunction has the form: 

*. = (5II^> (29) 

i<j 

consisting of a symmetrized product of pair correlation operators Fij operating 
on the Fermi gas wavefunction $. In symmetric nuclear matter, the function 
F^ includes eight terms: 

F i3 = E FiriMv (3°) 

P=l,8 

representing central, spin-spin, tensor and spin-orbit correlations with and 
without isospin factors. In pure neutron matter, the F^ reduce to a sum of four 
terms with only odd p <7. The correlation operators F^ are determined from 
Euler-Lagrange equations [67] that minimize the two-body cluster contribution 
of an interaction (V — A), where: 

V i3 = E c?VV{r l3 )0% (31) 

p=l,14 

K = E nm)oi. (32) 

p=l,8 

The variational parameters a p are meant to simulate the quenching of the 
spin-isospin interaction between particles i and j, due to flipping of the spin 
and/or isospin of particle i or j via interaction with other particles in matter. 
The NN interaction used in Refs. [12,13,66] is the recent parametrization of the 
Argonne group, the so-called Argonne Vis two-nucleon interaction discussed 
above. It has the form 

V^,ij= E V p {r l3 )O p l3 + V em . (33) 

p=l,18 

The electromagnetic part V em consists of Coulomb and magnetic interactions 
in the nn, np and pp pairs, and it is omitted from all nuclear matter studies. 
The strong interaction part of the potential includes fourteen isoscalar terms 



27 



with operators among these the central, spin-spin, spin-orbit and tensor op- 
erators. In addition, a phenomenological three-body force Vij k is included, 
represented by the Urbana models of V^ containing two isoscalar terms: 

Vw = V% + V t % . (34) 

The first term represents the Fujita-Miyazawa two-pion exchange interaction: 

Vijk = ^cyc Mir [{n ■ Tj, Ti ■ Tfc} {Xij, X ik } + \[t { - Tj,Ti • T k ] [X ij} X ik ])($5) 

Xij = SijT^rij) + Oi ■ <TjY„(rij), (36) 

with strength A 2n and where the matrices a and r are the Pauli matrices for 
spin and isospin, respectively. The functions T^{r\^) and Y^r^) describe the 
radial shapes of the one-pion exchange tensor and Yukawa potentials. These 
functions are calculated using the average value of the pion mass and include 
the short-range cutoffs used in the Argonne V\g NN interaction. The term 
denoted by V^ k is purely phenomenological, and has the form: 

V l % = U,Y.Tl{n j )Tl{r lk ). (37) 

eye 

This term is meant to represent the modification of NA- and A A-contributions 
in the two-body interaction by other particles in the medium, and also ac- 
counts for relativistic effects. The spin-isospin dependence of these effects is 
neglected. The two parameters A 2n and Uq are chosen to yield the observed 
energy of 3 H and the equilibrium density of nuclear matter, p — 0.16 fm~ 3 . 
Obviously, this fitting procedure will yield different parameters if another NN 
interaction is employed, e.g. the CD-Bonn interaction since, see discussion 
above, the various potentials yield slightly different binding energies for 3 H at 
the two-body level. The inclusion of many-body clusters is described in Ref. 
[62]. Relativistic boost corrections were also evaluated in Refs. [12,13,66] but 
the latter will be included in our discussion of relativistic effects in subsection 
2.2.3 below. 

In the remainder of this subsection we will henceforth discuss the consequences 
for the EoS from the above many-body corrections. 

In Fig. 13 we plot the results for PNM and SNM obtained with two-body 
interactions including only PPHH diagrams or higher-order diagrams stem- 
ming from the variational cluster approach. It is noteworthy to observe that 
in PNM the energy per particle up to 2-3 times no is rather similar for all cal- 
culations at the two-body level. To a certain extent this is expected since the 
strong nuclear tensor force contribution from the T = channel is absent. This 
means in turn that more complicated many-body terms are not so important 
in PNM. We will see this also in connection with the three-hole line discussion 
below. For symmetric nuclear matter the situation is however different due to 
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the strong tensor force in the np channel, leading to larger higher-order cor- 
rections. This is clearly seen in the lower panel of Fig. 13. In e.g. the PPHH 
calculation, where the correlations are due to hole-hole and particle-particle 
ladders, the tensor force plays the main role. Again this is similar to the situa- 
tion for the effective interaction in finite nuclei [42]. There, to e.g. second order 
in the interaction, diagrams with hole-hole and particle-particle intermediate 
states tend to be bigger than screening diagrams. In general however, none 
of the calculations at the two-body level reproduce properly the saturation 
properties of SNM. 

In Fig. 14 we have then included results from calculations with effective three- 
body terms and the phenomenological three-body forces. The results are com- 
pared with those from LOB with the Argonne Vis interaction. With the in- 
clusion of the phenomenological three-body force described in Eq. (34), which 
was fitted to reproduce the binding energy of the triton and the alpha particle, 
a clear change is seen both in PNM and SNM. The energy per particle gets 
more repulsive at higher densities. Three hole-line diagrams however, with the 
continuous choice are rather close to LOB with the continuous choice in SNM, 
while they are almost negligible in PNM, in line with our observation above 
about the tensor force component in the T = channel. 

Thus, as a summary, the inclusion of phenomenological three-body forces in 
non-relativistic calculations are needed in order to improve the saturation 
properties of the microscopically calculated EoS, whereas the inclusion of ad- 
ditional many-body effects, either three-hole line diagrams or through VCS 
calculations yield results at low densities similar to those from LOB. How- 
ever, even with the three-body force, one is not able to reproduce properly 
the binding energy. In the next subsection we will thence discuss further cor- 
rections stemming from relativistic effects. 



2.2.3 Relativistic effects 

The properties of neutron stars depend on the equation of state at densities up 
to an order of magnitude higher than those observed in ordinary nuclei. At such 
densities, relativistic effects certainly prevail. Among relativistic approaches 
to the nuclear many-body problem, the so-called Dirac-Hartree and Dirac- 
Hartree-Fock approaches have received much interest (Serot and Walecka [30]; 
Serot 1992 [68] Horowitz and Serot [69]). One of the early successes of these ap- 
proaches was the quantitative reproduction of spin observables, which are only 
poorly described by the non-relativistic theory. Important to these methods 
was the introduction of a strongly attractive scalar component and a repulsive 
vector component in the nucleon self-energy (Brockmann [70] and Serot and 
Walecka [30]). Inspired by the successes of the Dirac-Hartree- Fock method, a 
relativistic extension of Brueckner theory was proposed by Celenza and Shakin 
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Fig. 13. Upper panel: Energy per particle in PNM for various many-body approaches 
with two-body interactions only, i.e. PPHH results with the CD-Bonn interaction 
(stable results were obtained up to densities 0.4 fm -3 only), LOB with the Argonne 
Vis interaction and variational cluster (VCS) results for the Argonne Vis interaction 
model. Lower panel: The corresponding results for SNM. 

[71], known as the Dirac-Brueckner theory. One of the appealing features of 
the Dirac-Brueckner approach is the self-consistent determination of the rela- 
tivistic sp energies and wave functions. The Dirac-Brueckner approach differs 
from the Dirac-Hartree-Fock one in the sense that in the former one starts 
from the free NN interaction which is only constrained by a fit to the NN 
data, whereas the Dirac-Hartree-Fock method pursues a line where the pa- 
rameters of the theory are determined so as to reproduce the bulk properties 
of nuclear matter. It ought, however, to be stressed that the Dirac-Brueckner 
approach, which starts from NN interactions based on meson exchange, is 
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Fig. 14. Upper panel: Energy per particle in PNM for variational calculations 
with the Argonne Vis interaction without (VCS) and with three-body interactions 
(V18+UIX). The result from the three-hole line expansion of Baldo and co-workers 
employing the V14 interaction is included as well. Lower panel : The corresponding 
results for SNM. 

a non-renormalizable theory where the short-range part of the potential de- 
pends on additional parameters like vertex cut-offs, clearly minimizing the 
sensitivity of calculated results to short- distance inputs, see e.g. Refs. [71- 
73]. The description presented here for the Dirac-Brueckner approach follows 
closely that of Brockmann and Machleidt [72] . We will thus use the meson- 
exchange models of the Bonn group, denned in table A. 2 of Ref. [39]. There 
the three-dimensional reduction of the Bethe-Salpeter equation as given by 
the Thompson equation is used to solve the equation for the scattering matrix 
[74]. Hence, including the necessary medium effects like the Pauli operators 
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discussed in the previous subsection and the starting energy, we shall rewrite 
Eq. (23) departing from the Thompson equation. Then, in a self-consistent 
way, we determine the above-mentioned scalar and vector components which 
define the nucleon self-energy Note that negative energy solutions are not in- 
cluded. An account of these can be found in the recent work of de Jong and 
Lenske [75]. In order to introduce the relativistic nomenclature, we consider 
first the Dirac equation for a free nucleon, i.e. 

(i /) — m)ip(x) = 0, 

where m is the free nucleon mass and ip(x) is the nucleon field operator (x is a 
four-point) which is conventionally expanded in terms of plane wave states and 
the Dirac spinors u(p, s), and v (p, s), where p = (p°, p) is a four momentum^ 
and s is the spin projection. 

The positive energy Dirac spinors are (with uu = 1) 



u(p, s) 



1 E(p) + m 
2m 



' x, ^ 



<T p 

\ E(p)+m,Xs J 



(3* 



where Xs is the Pauli spinor and E(p) 



in- 



+ |p| 2 . To account for medium 



modifications to the free Dirac equation, we introduce the notion of the self- 
energy S(p). As we assume parity to be a good quantum number, the self- 
energy of a nucleon can be formally written as 



S(p) = S 5 (p)-7oS°(p)+7pS y (p). 

The momentum dependence of S° and S^ is rather weak (Serot and Walecka 
1986). Moreover, TX « 1, such that the features of the Dirac-Brueckner- 
Hartree-Fock procedure can be discussed within the framework of the phe- 
nomenological Dirac-Hartree ansatz, i.e. we approximate 



S « S 5 - 7o S° = U s 



U v , 



where Us is an attractive scalar field, and Uy is the time-like component of a 
repulsive vector field. The finite self-energy modifies the free Dirac spinors of 



Further notation is as given in Itzykson and Zuber [76]. 
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Eq. (38) as 



u(p, s) 



E(p) + m 
2rh 



Xs 



<7P v 

\ E(p)+m Xs ) 



where we let the terms with tilde represent the medium modified quantities. 
Here we have defined 



m = m + Us, 



and 



Ei = E( Pl 



m\ + pi 



(39) 



As in the previous subsection, the subscripts i, and h below, represent the 
quantum numbers of the single-particle states, such as isospin projections Tj 
and Th, momenta ki and kh, etc. 



The sp energy is 



Si = Ei + U 



v- 



and the sp potential is given by the G-matrix as 



U; 



Y, ^(ih\G(E = E i + E h )\ih) 



h<k F EiE h 



AS ' 



(40) 



(41) 



or, if we wish to express it in terms of the constants Us and Uy, we have 
rhi 



U; 



-^U S * + Uy. 



(42) 



In Eq. (41), we have introduced the relativistic G-matrix. If the two interact- 
ing particles, with isospins T\ and r 2 , give a total isospin projection T z , the 
relativistic G-matrix in a partial wave representation is given by 



Gf*{kk'KE) = V l f>{kk') + Y J f 



d 3 q 



(43) 



raTzn ^rh 2 Q^( q ,K) _ d «T f{qk , K ^ 



XV "" Z{kq) E!EHE-El-El) 
where the relativistic starting energy is defined according to Eq. (39) as 



E = E(Jk 2 + K 2 /4, n) + E( Jk* + K 2 /4, r 2 ). 
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and 

Equations (40)- (43) are solved self-consistently, starting with adequate values 
for the scalar and vector components Us and Uy. This iterative scheme is 
continued until these parameters show little variation. Finally, the relativistic 
version of Eq. (26) reads 

£ /A = 1 m^ + (44) 

71 h<k F &h 

±- £ ^^(hti\G(E = E h + E h >)\hh') AS -m. 



An alternative approach which we will also discuss is the inclusion of rela- 
tivistic boost corrections to non-relativistic NN interactions in the VCS cal- 
culations of Akmal et al. [13,66]. In all analyses, the NN scattering data are 
reduced to the center-of-mass frame and fitted using phase shifts calculated 
from the nucleon-nucleon interaction, V, in that frame. The interaction ob- 
tained by this procedure describes the NN interaction in the cm. frame, in 
which the total momentum P^ = p« + Pj, is zero. In general, the interaction 
between particles depends upon their total momentum, and can be written as 

v(P ij ) = v ij + 8v(P ij ), (45) 

where i>„ is the interaction for P^- = 0, and 5v(Pij) is the boost interaction 
[77] which is zero when P^- = 0. 

Following the work of Krajcik and Foldy [78], Friar [79] obtained the following 
equation relating the boost interaction of order P 2 to the interaction in the 
center of mass frame: 

5v(P) = --^v + -L[P • r P • V, v] + -^[to - ffj) x P • V,v]. (46) 
8m 2 8m 2 8m 2 

The general validity of this equation in relativistic mechanics and field theory 
was recently discussed [77] . Incorporating the boost into the interaction yields 
a non-relativistic Hamiltonian of the form: 

h nr = Et^ + £(% + MPy)) + E V3* + • • • , (47) 

where the ellipsis denotes the three-body boost, and four and higher body 
interactions. This H^ R contains all terms quadratic in the particle velocities, 
and is therefore suitable for complete studies in the non-relativistic limit. 
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The authors of Refs. [64,80] find that the contribution of the two-body boost 
interaction to the energy is repulsive, with a magnitude which is 37% of the 
VM. contribution. The boost interaction thus accounts for a significant part of 
the V^ k in Hamiltonians which fit nuclear energies neglecting Sv. 

In this work we will follow Refs. [12,13,66] and only keep the terms of the boost 
interaction associated with the static part of V, and neglect the last term 
in Eq. (46). That term is responsible for Thomas precession and quantum 
contributions that are negligibly small here [81]. The correction Sv is then 
given by 

Sv(P) = ~^v s + -Lp • r P • W. (48) 

The two terms are due to the relativistic energy expression and Lorentz con- 
traction, and are denoted Sv RE and Sv LC , respectively. The three-nucleon in- 
teraction used in the H* NR of Eq. (47) is denoted by V*j k . Its parameters are 
obtained by fitting the binding energies of 3 H and 4 He, and the equilibrium 
density of SNM, including Sv. The strength of V R k is 0.63 times that of V R k 
in UIX, while that of V t jl is unchanged. The resulting model of Vij k is called 
UIX*. 

In Fig. 15 we plot the non-relativistic BHF and DBHF energies per particle 
for both SNM and PNM. These results were discussed in Ref. [45] and were 
obtained with an older version of the Bonn interactions, namely the so-called 
Bonn A potential. The relativistic corrections include in an effective way a 
certain class of three-body corrections, the so-called Z-graphs [82], although 
the calculations are done only at the two-body level. The non-relativistic re- 
sults are similar to those obtained with the CD-Bonn interaction discussed in 
the previous subsection. One notices also that the DBHF results yield satu- 
ration properties close to the empirical data, a binding energy of —15.1 MeV 
and saturation density of 0.19 fm~ 3 , while the non-relativistic results saturate 
at 0.37 fm~ 3 with a binding energy of —22.4 MeV. For comparison we have 
also included the results of Akmal et al. [12,13] with three-body diagrams and 
boost corrections. Again, one notices that at densities below tiq all many-body 
schemes tend to give similar results, although in this case, the relativistic re- 
sults are in better agreement with the results of the full calculation of Akmal 
et al. [12,13]. For completeness, we display the energy per particle for /3-stable 
matter obtained with the DBHF and the results from the calculations of Ak- 
mal et al. [13,66] with boost corrections and three-body forces in Fig. 16. The 
corresponding proton fractions are also plotted. The dip in the proton fraction 
of Akmal et al. at ~ 0.2 fm~ 3 is due to the formation of a pion condensate. 
For higher densities we note that BHF and the calculations of Akmal et al. 
yield rather similar proton fractions. The direct Urea starts only at densities 
~ 1 fm~ 3 , while the DBHF calculations allow the direct Urea at much lower 
densities. 
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Fig. 15. Upper panel: DBHF and BHF energy per nucleon for PNM. The results with 
boost corrections of Akmal et al. [13] with (Fig + (5-u+UIX*) and without (Fig + 5v) 
three-body forces. Lower panel: the corresponding results for SNM. 

2.3 A causal parametrization of the nuclear matter EoS 



Since three-body forces are expected to be important (it suffices to mention 
studies of the Triton [44]), we will in our discussions of the mixed phase in the 
next section and in connection with the structure of a neutron star, employ 
the recent EoS of Akmal et al. [13]. A non-relativistic EoS (although with the 
inclusion of boost corrections) is preferred here. The EoS for nuclear matter is 
thus known to some accuracy for densities up to a few times nuclear saturation 
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Fig. 16. Upper panel: Energy per nucleon in (3 — stable matter for the DBHF ap- 
proach and the results of Akmal et al. with boost corrections and three-body forces 
(Vis + Sv+VIX*). Lower panel: the corresponding proton fraction x p . For the DBHF 
calculation both electrons and muons are included. 



density, n = 0.16 fm~ 3 . Detailed knowledge of the EoS is crucial for the 
existence of, e.g. pion condensates [12,13] or the delicate structures in the 
inner crust of neutron stars [26]. However, for the gross properties and our 
discussion of properties of neutron stars we will adopt a simple form for the 
binding energy per nucleon in nuclear matter consisting of a compressional 
term and a symmetry term 



£ = E comp (n) + S(n)(l - 2x P Y 



SqU 



u 



s 



l + 5u 



+ Soii^l - 2x P 



(49) 
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Here u = n/no is the ratio of the baryon density to nuclear saturation density 
and we have defined the proton fraction x v = n p /n. The compressional term 
is in (49) parametrized by a simple form which reproduces the saturation 
density, binding energy and compressibility. The binding energy per nucleon 
at saturation density excluding Coulomb energies is So = — 15.8MeV and the 
parameter S = 0.2 is determined by fitting the energy per nucleon at high 
density to the EoS of Akmal et al. [13] with three-body forces and boost 
corrections, but taking the corrected values from Table 6 of Ref. [13]. The 
reason behind the choice of 5 = 0.2 will explained below in connection with 
the discussion of Figs. 18 and 19. Further discussions will be presented in 
section 5. The corresponding compressibility is K = 18£o/(l + $) — 200MeV 
in agreement with the experimental value. For the symmetry term we obtain 
So = 32 MeV and 7 = 0.6 for the best fit. The quality of this simple functional 
exhibits a x 2 P er datum close to 1 and is compared with the results of Akmal 
et al. [13] for both PNM and SNM in Fig. 17. As can been from this figure, 
the agreement is rather good except at the very high densities where the 
EoS of Akmal et al. [13] becomes superluminal and therefore anyway must 
be wrong. A much more sophisticated fit which reproduces the data in terms 
of Skyrme functional approach is given by the Akmal et al. [13]. However, it 
is amazing that such a simple quadratic formula fits so well the data coming 
from a microscopic calculation. In view of the uncertainties which pertain 
to the EoS at higher densities, we feel that our parametrization is within 
present error margins. The agreement between the microscopic calculation of 
Akmal et al. [13] and the simple parametrization of Eq. (49) may imply that 
the essential many-body physics close to the equilibrium density arises from 
two-body and three-body terms to S only. The reason beingP] that three-body 
terms are proportional with n 3 while the two-body terms are proportional with 
n 2 . With three-body terms we obviously intend both effective interactions and 
contributions from real three-body forces. The evaluation of the latter is, as 
discussed above, still an unsettled problem. 

If we now restrict the attention to matter with electrons only, one can easily 
obtain an analytic equation for the proton fraction through the asymmetry 
parameter x. Recalling the equilibrium conditions for /^-stable matter of Eqs. 
(6)-(7) and using the definitions of the chemical potentials for particle species 
i of Eq. (3) one finds that 

1 de 
* = ndx~ p > (50) 

and from the latter it is rather easy to show that the asymmetry parameter x 



This argument is for the energy density, i.e. e = En. 
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is given by (assuming ultra-relativistic electrons) 



nx r 



(4,W(1 - 2x p )f 
3vr 2 



(51) 



Defining 



2 (ASqv? 
n 2 n 



(52) 
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Eq. (51) reduces to 

3x 3 + ax — a = , (53) 

where x = 1 — 2x p . Since we will always look at solutions for densities greater 
than zero, the cubic equation for x has actually an analytical solution which 
is real and given by 

i 
with tanip = (tan | J 3 and tan0 = —2y/a/3. Note well that x depends on the 
total baryon density n only. This means in turn that our parametrization of 
the EoS can now be rewritten for /3-stable matter as 

and is an analytical function of density only. 

The quality of our approximation to the EoS of Akmal et al. [13] for other 
observables than the energy per particle is shown in Fig. 18 for the proton 
fractions derived from the simple expression in Eq. (54). In the same figure we 
display also the resulting energy per nucleon in /3-stable matter and compare 
it with the results of Akmal et al. [13] for various values of 5. Note well that 
the proton fraction does not depend on the value of 5, see Eq. (51). As can 
be seen from this picture, the EoS with 5 = yields the stiff est EoS, and as 
a consequence it results in a superluminal behavior at densities greater than 
n » 1.0 fm~ 3 . This is seen in Fig. 19 where we plot the sound speed (v s /c) 2 
for various S values and that resulting from the microscopic calculation of 
Akmal et al. [13] The form of (49), with the inclusion of the parameter 5, 
provides therefore a smooth extrapolation from small and large densities with 
the correct behavior in both limits, i.e. the binding energy per nucleon E/A = 
S is linear in number density. In the dilute limit this is the Lenz (optical) 
potential. At high densities the linearity is required by the condition that the 
sound speed c 2 s = dP/de does not exceed the speed of light. This justifies the 
introduction of the parameter 5 in our parametrization and explains also our 
deviation from the results of Akmal et al. at densities greater than 0.6 ~ 0.7 
fm~ 3 , see Figs. 17 and 18. For S = 0.1 the EoS becomes superluminal at 
densities of the order of 6 fm -3 . From the definition of (v s /c) 2 in Eq. (18) and 
the EoS of Eq. (49) that for 



So 




8>J—~0.13, (56) 

V m n 

the EoS remains causal for all densities. The EoS of Akmal et al. becomes 
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Fig. 18. Upper panel: Energy per nucleon in /3-stable matter for the parametrized 
EoS of Eq. (49) for 6 = 0, 0.13, 0.2, 0.3 and the results of Akmal et al. [13] with boost 
corrections and three-body forces (Vis + Sv+VTK*). Lower panel: the corresponding 
proton fraction x p . 

superluminous at n ps 1.1 fm -3 . With this caveat we have an EoS that re- 
produces the data of Akmal et al. at densities up to 0.6 ~ 0.7 fm -3 and has 
the right causal behavior at higher densities. Furthermore, the differences at 
higher densities will also not be of importance in our analysis of the dynamics 
and structure of neutrons stars, since the mixed baryon-quark phase, with 
realistic values for the bag parameter and the coupling constant a s starts at 
densities around 0.5 — 0.8 fm~ 3 . 

Finally, in Fig. 19 we have also plotted the sound speed following the ap- 
proach of Baym and Kalogera [83], where the sound speed is allowed to jump 
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Fig. 19. (v s /c) 2 for 5 = 0, 0.1, 0.2, 0.3, the results of Akmal et al. [13], and for Baym 
and Kalogera's [83] patched EoS which shows a discontinuous (v s /c) 2 . 

discontinuously at a chosen density in order to keep the EoS causal. With this 
prescription, Baym and Kalogera were also able to obtain an upper bound for 
neutron star masses of 2.9M . The approach of Baym and Kalogera differs 
substantially from ours since their EoS is discontinuously stiffened by taking 
v s — c at densities above a certain value n c which, however, is lower than 
n s = 5hq where their nuclear EoS becomes superluminous. This stiffens the 
nuclear EoS for densities n c < n < n s but softens it at higher densities. Their 
resulting maximum masses lie in the range 2.2M Q < M < 2.9M Q . Our ap- 
proach incorporates causality by reducing the sound speed smoothly towards 
the speed of light at high densities. Therefore our maximum mass never ex- 
ceeds that the of nuclear EoS of Akmal et al. [13]. In fact one may argue that 
at very high densities particles become relativistic and the sound speed should 
be even lower, v 2 ~ c 2 /3, and therefore the softening we get from incorporating 
causality is even on the low side. 



2.4 Hyperonic matter 



At nuclear matter density the electron chemical potential is ~ 110 MeV, see 
e.g. Fig. 6. Once the rest mass of the muon is exceeded, it becomes energetically 
favorable for an electron at the top of the e~ Fermi surface to decay into a 
H~ . We then develope a Fermi sea of degenerate negative muons, see again 
Fig. 6. In a similar way, as soon as the chemical potential of the neutron 
becomes sufficiently large, energetic neutrons can decay via weak strangeness 
non- conserving interactions into A hyperons leading to a A Fermi sea with 
^K — ^n- However, if we neglect interactions, or assume that their effects are 
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Fig. 20. The neutron and electron chemical potentials in beta stable matter ac- 
cording to models V±8 + fo+UIX* (full line) and Vis + $v (dashed line). Threshold 
densities for the appearance of non-interacting hyperons are marked by horizontal 
line segments. Taken from Ref. [13]. 

small, one would expect the S~ to appear via 



n 



v, 



(57) 



at lower densities than the A, even though S~ is more massive the reason 
being that the above process removes both an energetic neutron and an en- 
ergetic electron, whereas the decay to a A, being neutral, removes only an 
energetic neutron. Stated differently, the negatively charged hyperons appear 
in the ground state of matter when their masses equal \x e + //„, while the neu- 
tral hyperon A appears when \x n equals its mass. Since the electron chemical 
potential in matter is larger than the mass difference ms- -wa — 81.76 MeV, 
the YT will appear at lower densities than the A. We show this in Fig. 20 
where we plot the chemical potentials for electrons and neutrons in /3-stable 
matter. The threshold densities for E~, A and the isobar A - are indicated by 
the horizontal lines. 
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Since this work has an emphasis on many-body approaches, we will try to 
delineate here as well how to obtain properties of hyperons in dense nuclear 
matter within the framework of microscopic theories. The main problem we 
have to face in our case is that the hyperon-nucleon interaction and especially 
the hyperon-hyperon interaction are less constrained by the data as is the 
case in the nucleonic sector. For a recent version of the Nijmegen model for 
the baryon-baryon interaction and applications to nucler matter, see Stoks et 
al. in Refs. [18,19]. 

In spite of these problems, we will attempt at describing how to approach hy- 
peronic matter starting with realistic hyperon-nucleon and hyperon-hyperon 
interactions. We will start with the simplest possible many-body scheme again, 
namely lowest order Brueckener theory (LOB). The nuclear matter YN G- 
matrix is solved in momentum space and the two-particle YN states are de- 
fined in terms of relative and the center-of-mass momenta, k and K, given 
by 



M N k Y - M y k 



N 



k 

M N + M Y 
K = kjv + k y , 

where k^ and ky are the nucleon and hyperon momenta, respectively. Using 
an angle-averaged Pauli operator, we perform a partial wave decomposition 
of the Bethe-Goldstone equation, which, in terms of the quantum numbers of 
the relative and center-of-mass motion (RCM), is written as 

((Y"N)k"l"KS"\G aT *{uj) \{YN)klKS) = 
{{Y"N)k"l"KS"\ V aT * \{YN)klKS) 

+E E / k ' 2dk> ( ( Y " N ) k"l"KS" | V aT * | (Y'N) k'l'KS') 

US' Y'=AE J 

Q YN (k',K) 



X 



~ kV(M N + M ,) 

^ * M, +M Y 



2{M n +M y ,) 2M n M y , ±y± Y' 



x((Y'N)k'l'KS'\G aT *(uj)\(YN)klKS} , (58) 

where Qyn is the hyperon-nucleon nuclear Pauli operator, V is the YN po- 
tential and uj is the starting energy which corresponds to the sum of non- 
relativistic single-particle energies of the interacting nucleon and hyperon. 
Note that kinetic energies are used in the intermediate Y'N states and My — 
My accounts for the mass difference of the initial and intermediate hyperon. 
The label a represents in this the angular momentum J and the orbital an- 
gular momentum L of the center-of-mass (cm.) motion. The variables k, k', 
k" and /, I', I" denote relative momenta and angular momenta, respectively, 
while K is the momentum of the center of mass motion. 
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In the Brueckner-Hartree-Fock approach the hyperon single-particle potential 
U Y is obtained self-consistently by the following sum of diagonal G-matrix 
elements: 

Uy (ky) = J d 3 k N (Yk Y , Nk N \ G(e N (k N )+e Y (k Y )) \Yk Y , Nk N ) , (59) 

where e N (y)(kN(Y)) = k 2 NiY J \2M N ^ Y y) + U N ( Y )(k N ( Y )) is the single-particle 
energy of the nucleon (hyperon). The superscript iV refers to the fact that 
we have only nucleon holes. Using the partial wave decomposition of the G- 
matrix, the single-particle potential U Y can be rewritten as 

^^^TTn E (2J + 1)(2T+1) 

x f k 2 dkf(k,k Y ) (YNklS\G aT *(cu)\YNklS) , (60) 

o 

where an average over the hyperon spin and isospin (t Y ) has been performed 
and the weak center-of-mass dependence of the G-matrix has been neglected. 
In Eq. (60), k is the relative momentum of the YN pair, £y = Mn/M y , k max 
is given by 

kmax = '-, : I (61) 

and the weight function f(k, ky) by 



f(k,h 



1 for k < ^gg^ , 

for \£ Y k Y ~ (1 + fr)k\ > k F , (62) 



k^[( Y k Y -(iHv)kV otherwige 
4( Y (i+(Y)k Y k o^neiwise 



The binding energy of the hyperon to its single-particle energy is then 

B Y {k Y )^e{k Y ) = ^- + U^{k Y ) (63) 

However, we see from Eq. (59) that we have only nucleons as intermediate hole 
states. This approach is viable only if we are interested in studying properties 
of hyperons on top of /3-stable matter. Such studies were recently performed 
by Schulze et al. [14]. These calculations could be viewed as way of establish- 
ing the possibility of hyperon formation in matter. The results for /5-stable 
matter with hyperon formation from Ref. [14] is shown in Fig. 21 employ- 
ing various models for the EoS. The soft-core nucleon-hyperon interaction of 
the Nijmegen group is included here [85] in order to evaluate the hyperon 
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Fig. 21. The onset of £~ and A formation employing the Paris interaction in an 
BHF calculation, the Paris interaction added three-body forces (TBF), the Argonne 
V14 plus TBF from Ref. [84], and the DBHF calculations discussed in subsection 
2.2.3. The short dashes indicate the occurence of chemical equilibrium for T,~ and 
the long dashes the corresponding one for A. 

single-particle energy. Since the EoS which includes three-body interactions 
or relativistic effects is stiffer than the non-relativistic EoS from BHF calcu- 
lations (the results labelled Paris), one expects the threshold for formation to 
be lower. We note from the above figure that the relativistic and three-body 
results seem to agree fairly well and yield an onset of hyperon formation at 
densities ~ 2 — 3n . 

However, with a finite hyperon concentration, this entails also a summation 
over hyperon holes in Eq. (59). Eq. (63) gets an additional term Uyiky) where 
the latter is obtained from an interaction between hyperons only. Furthermore, 
in order to obtain the correct value for the Fermi surface in (3-stable matter 
with hyperons, one needs to solve e.g. Eqs. (8) and (9) as well. Results of such 
calculations are in progress [86] with realistic hyperon-hyperon interactions. 

In this work, we will make life somewhat easier by just including quark degrees 
of freedom through a simplified model, namely the bag model, in order to 
account for degrees of freedom beyond the nucleonic ones. This is discussed in 
subsection 2.8. 
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2.5 Kaon condensation 



Kaon condensation in dense matter was suggested by Kaplan and Nelson [11], 
and has been discussed in many recent publications [87,88]. Due to the attrac- 
tion between K~ and nucleons its energy decreases with increasing density, 
and eventually if it drops below the electron chemical potential in neutron star 
matter in /^-equilibrium, a Bose condensate of K~ will appear. It is found that 
K~'s condense at densities above ~ 3 — 4po, where po — 0.16 fm~ 3 is normal 
nuclear matter density. This is to be compared to the central density of ~ 4p 
for a neutron star of mass 1.4M according to the estimates of Wiringa, Fiks 
and Fabrocini [84] using realistic models of nuclear forces. 

In neutron matter at low densities, when the interparticle spacing is much 
larger than the range of the interaction, tq 3> R, the kaon interacts strongly 
many times with the same nucleon before it encounters and interacts with 
another nucleon. Thus one can use the scattering length as the "effective" 
kaon-nucleon interaction, a^-N — — 0.41fm. The kaon energy deviates from 
its rest mass by the Lenz potential 

2tt . . 

ULenz = m K H a K - N p, (64) 

m R 

which is the optical potential obtained in the impulse approximation. If hadron 
masses furthermore decrease with density the condensation will occur at lower 
densities [87]. 

At high densities when the interparticle spacing is much less than the range 
of the interaction, tq <C R, the kaon will interact with many nucleons on a 
distance scale much less the range of the interaction. The kaon thus experiences 
the field from many nucleons and the kaon energy deviates from its rest mass 
by the Hartree potential: 

U Hartree = Vd K + P V K - N (r)d 3 r , (65) 

As shown in Ref . [89] , the Hartree potential is considerably less attractive than 
the Lenz potential. Already at rather low densities, when the interparticle dis- 
tance is comparable to the range of the KN interaction, the kaon-nucleon 
and nucleon-nucleon correlations conspire to reduce the K~N attraction sig- 
nificantly [90]. This is also evident from Fig. 2.5 where the transition from 
the low density Lenz potential to the high density Hartree potential is calcu- 
lated by solving the Klein-Gordon equation for kaons in neutron matter in the 
Wigner-Seitz cell approximation. Results are for square well i^~iV-potentials 
of various ranges R. 

Kaon-nucleon correlations reduce the K~N interaction significantly when its 
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Fig. 22. Kaon energy as function of neutron density. Including nuclear correlations 
in the Wigner-Seitz cell approzimation is shown by full curves for varios ranges of 
the K~n potentials R = 0.4 fm, R = 0.7 fm and R = 1.0 fm. At low densities they 
approach the Lenz result (Eq. (64), dotted curve) and at high densities they ap- 
proach the Hartree result (Eq.(65), dashed curves). The electron chemical potential 
H e of our EoS (Eq. 49) with 5 = 0.2 is shown with (lowest dotted curve) and without 
a transition to a mixed phase of quark matter for a Bag constant of B 1 ' 4 = 100 
MeVfrn" 3 . 

range is comparable to or larger than the nucleon-nucleon interparticle spac- 
ing. The transition from the Lenz potential at low densities to the Hartree 
potential at high densities occurs already well below nuclear matter densities. 
For the measured K~n scattering lengths and reasonable ranges of interactions 
the attraction is reduced by about a factor of 2-3 in cores of neutron stars. Rel- 
ativistic effects further reduce the attraction at high densities. Consequently, 
a kaon condensate is less likely in neutron stars due to nuclear correlations. 
However, if kaon masses drop with densities [87] condensation will set in at 
lower densities. 

If the kaon condensate occurs a mixed phase of kaon condensates and ordinary 
nuclear matter may coexist in a mixed phase [91] depending on surface and 
Coulomb energies involved. The structures would be much like the quark and 
nuclear matter mixed phases described above. 
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2.6 Pion condensation 



Pion condensation is like kaon condensation possible in dense neutron star 
matter. For an in depth survey see e.g. Refs. [35,92] and references therein. 
If we first neglect the effect of strong correlations of pions with the matter in 
modifying the pion self-energy, one finds it is favorable for a neutron on the 
top of the Fermi sea to turn into a proton and a 7r~ when 

«n = Up + He > "V, (66) 

where m n = 139.6 MeV is the n~ rest mass. As discussed in the previous sub- 
section, at nuclear matter saturation density the electron chemical potential is 
~ 100 MeV and one might therefore expect the appeareance of tt~ at a slightly 
higher density. One can however not neglect the interaction of the pion with 
the background matter. Such interactions can enhance the pion self-energy 
and thereby the pion threshold density, and depending on the chosen param- 
eters , see again Ref. [35], the critical density for pion condensation may vary 
from no to 4uq. These matters are however not yet settled in a satisfying way, 
and models with strong nucleon-nucleon correlations tend to suppress both the 
nNN and irAN interaction vertices so that a pion condensation in neutron 
star matter does not occur. However, in addition to a charged pion conden- 
sate, one may also form a it condensate through the reaction n — ► n + 7r° 
if the 7T° effective mass in the medium is zero. The recent analysis, based on 
the Vis interaction model [37] of Akmal et al. [12,13,66] suggests such a pion 
condensate. The effects were partly discussed in section 2.2 and the impact on 
the proton fraction was shown in Fig. 16. The 7r° condensation of Akmal et 
al. [13] for pure neutron matter appears at a density of ~ 0.2 fm~ 3 when the 
three-body interaction is included, whereas without V^ it appears at much 
higher densities, i.e. ~ 0.5 fm~ 3 . Although it is a robust mechanism in the 
variational calculation of Ref. [13], the conclusion relies on both the model 
of the NN and three-body interactions adopted in the calculations. As noted 
in Ref. [40] for pure neutron matter, the Vis interaction resulted in a slightly 
different energy per particle at densities greater than 0.4 fm -3 when compared 
with the CD-Bonn and Nijmegen interactions. This topic was also discussed 
in in section 2.2 and shown in Fig. 3. Thus, before a firm conclusion can be 
reached about 7r° condensation, it is our belief that it should also be obtained 
at the two-body level with the other phase-shift equivalent NN interactions. 
That would lend strong support to the conclusions reached in Ref. [13]. The 
inclusion of three-body interactions introduces a further model dependence. 

Due to these uncertainties, we will refrain in this work from presenting a 
thourough discussion of pion condensation. Rather, we will take the liberty to 
refer to, e.g. Refs. [13,28,35,92]. 
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2.7 Superfluidity in baryonic matter 



The presence of neutron superfluidity in the crust and the inner part of neutron 
stars are considered well established in the physics of these compact stellar ob- 
jects. In the low density outer part of a neutron star, the neutron superfluidity 
is expected mainly in the attractive 1 S'o channel. At higher density, the nuclei 
in the crust dissolve, and one expects a region consisting of a quantum liquid 
of neutrons and protons in beta equilibrium. The proton contaminant should 
be superfluid in the 1 S channel, while neutron superfluidity is expected to 
occur mainly in the coupled 3 P2- 3 -^2 two-neutron channel. In the core of the 
star any superfluid phase should finally disappear. 

The presence of two different superfluid regimes is suggested by the known 
trend of the nucleon-nucleon (NN) phase shifts in each scattering channel. 
In both the 1 S'o and 3 P2- 3 -^2 channels the phase shifts indicate that the NN 
interaction is attractive. In particular for the 1 5'o channel, the occurrence of 
the well known virtual state in the neutron-neutron channel strongly suggests 
the possibility of a pairing condensate at low density, while for the 3 P2- 3 F 2 
channel the interaction becomes strongly attractive only at higher energy, 
which therefore suggests a possible pairing condensate in this channel at higher 
densities. In recent years the BCS gap equation has actually been solved with 
realistic interactions, and the results confirm these expectations. 

The 1 S'o neutron superfluid is relevant for phenomena that can occur in the 
inner crust of neutron stars, like the formation of glitches, which may to be 
related to vortex pinning of the superfluid phase in the solid crust [93] . The 
results of different groups are in close agreement on the 1 So pairing gap values 
and on its density dependence, which shows a peak value of about 3 MeV 
at a Fermi momentum close to kp ~ 0.8 fm _1 [94-97]. All these calculations 
adopt the bare NN interaction as the pairing force, and it has been pointed out 
that the screening by the medium of the interaction could strongly reduce the 
pairing strength in this channel [97-99] . However, the issue of the many-body 
calculation of the pairing effective interaction is a complex one and still far 
from a satisfactory solution. 

The precise knowledge of the 3 P2- 3 F2 pairing gap is of paramount relevance for, 
e.g. the cooling of neutron stars, and different values correspond to drastically 
different scenarios for the cooling process [28]. Generally, the gap suppresses 
the cooling by a factor ~ exp(— A/T), see e.g. Ref. [46], which is severe for 
temperatures well below the gap energy. Unfortunately, only few and partly 
contradictory calculations of the pairing gap exist in the literature, even at 
the level of the bare NN interaction [100-104] . However, when comparing the 
results, one should note that the NN interactions used in these calculations are 
not phase-shift equivalent, i.e. they do not predict exactly the same NN phase 
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shifts. Furthermore, for the interactions used in Refs. [100-103] the predicted 
phase shifts do not agree accurately with modern phase shift analyses, and the 
fit of the NN data has typically \ 2 /datum ps 3. As we discussed in subsection 
2.2, progress has been made not only in the accuracy and the consistency 
of the phase-shift analysis, but also in the fit of realistic NN interactions to 
these data. As a result, several new NN interactions have been constructed 
which fit the world data for pp and np scattering below 350 MeV with high 
precision. Potentials like the recent Argonne Vis [37], the CD-Bonn [38] or the 
new Nijmegen potentials [36] yield a x 2 /datum of about 1 and may be called 
phase-shift equivalent. In Table 1 we show the recent non-relativistic pairing 
gaps for the 3 P2- 3 -^2 partial waves, where effective nucleon masses from the 
lowest-order Brueckner-Hartree-Fock calculation of subsection 2.2 have been 
employed, see Ref. [105] for more details. These results are for pure neutron 
matter and we observe that up to kp ~ 2 fin -1 , the various potentials give more 
or less the same pairing gap. Above this Fermi momentum, which corresponds 
to a lab energy of ~ 350 MeV, the results start to differ. This is simply 
due to the fact that the potentials are basically fit to reproduce scattering 
data up to this lab energy. Beyond this energy, the potentials predict rather 
different phase shifts for the 3 P 2 - 3 -^2 partial waves, see e.g. Ref. [105]. Thus, 

Table 1 

Collection of 3 P2- 3 -^2 energy gaps (in MeV) for the various potentials discussed in 
subsection 2.2. BHF single-particle energies have been used. In case of no results, a 
vanishing gap was found. 



k F (fin -1 ) 


CD-Bonn 


Vis 


Nijm I 


Nijm II 


1.2 


0.04 


0.04 


0.04 


0.04 


1.4 


0.10 


0.10 


0.10 


0.10 


1.6 


0.18 


0.17 


0.18 


0.18 


1.8 


0.25 


0.23 


0.26 


0.26 


2.0 


0.29 


0.22 


0.34 


0.36 


2.2 


0.29 


0.16 


0.40 


0.47 


2.4 


0.27 


0.07 


0.46 


0.67 


2.6 


0.21 




0.47 


0.99 


2.8 


0.17 




0.49 


1.74 


3.0 


0.11 




0.43 


3.14 



before a precise calculation of 3 P2- 3 -^2 energy gaps can be made, one needs NN 
interactions that fit the scattering data up to lab energies of ~ 1 GeV. This 
means in turn that the interaction models have to account for the opening of 
inelasticities above 350 MeV due to the NA channel. 

The reader should however note that the above results are for pure neutron 
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matter. We end therefore this subsection with a discussion of the pairing gap 
for /3-stable matter of relevance for the neutron star cooling discussed in section 
5. We will also omit a discussion on neutron pairing gaps in the 1 Sq channel, 
since these appear at densities corresponding to the crust of the neutron star. 
The gap in the crustal material is unlikely to have any significant effect on 
cooling processes [26] , though it is expected to be important in the explanation 
of glitch phenomena. Therefore, the relevant pairing gaps for neutron star 
cooling should stem from the the proton contaminant in the 1 So channel, and 
superfluid neutrons yielding energy gaps in the coupled 3 P2- 3 F 2 two-neutron 
channel. If in addition one studies closely the phase shifts for various higher 
partial waves of the NN interaction, one notices that at the densities which 
will correspond to the core of the star, any superfluid phase should eventually 
disappear. This is due to the fact that an attractive NN interaction is needed 
in order to obtain a positive energy gap. 

Since the relevant total baryonic densities for these types of pairing will be 
higher than the saturation density of nuclear matter, we will account for rel- 
ativistic effects as well in the calculation of the pairing gaps. To do so, we 
resort to the Dirac-Brueckner-Hartree-Fock (DBHF) formalism discussed in 
subsection 2.2.3. 

As an example, consider the evaluation of the proton 1 So pairing gap using the 
DBHF approach. To evaluate the pairing gap we follow the scheme of Baldo et 
al. [94]. These authors introduced an effective interaction V^y- This effective 
interaction sums up all two-particle excitations above a cutoff momentum }cm, 
ku = 3 fm _1 in this work. It is defined according to 

Vk,k' = Vk,k' — 2^ VJ^fc"— — Vy t y, (67) 

k">k M tk " 



where the energy £ fc is given by 8k = y (£& — £>) + A^, e> being the single- 
particle energy at the Fermi surface, V^y is the free nucleon-nucleon potential 
in momentum space, defined by the three-momenta k, k'. The renormalized 
potential V^y and the free NN interaction V^y carry a factor rh 2 / E^Ey , due 
to the normalization chosen for the Dirac spinors in nuclear matter. These 
constants are also included in the evaluation of the G-matrix, as discussed in 
[42,72]. For the 1 S'o channel, the pairing gap A& is [94,106,107] 

A* 



k'<kM 



28y 



These equations are solved self-consistently in order to obtain the pairing gap 
A for protons and neutrons for different partial waves. 

In Fig. 23 we plot as function of the total baryonic density the pairing gap 



52 



for protons in the l So state, together with the results from the non-relativistic 
approach discussed in Refs. [103,109]. The results in the latter references were 
also obtained with the Bonn A potential of Ref. [39]. These results are all 
for matter in /^-equilibrium. In Fig. 24 we plot the corresponding relativistic 
results for the neutron energy gap in the 3 P2 channel. For the 1 D2 channel we 
found both the non-relativistic and the relativistic energy gaps to vanish. The 
non-relativistic results for the Bonn A potential are taken from Ref. [108]. 

As can be seen from Fig. 23, there are only small differences (except for higher 
densities) between the non-relativistic and relativistic proton gaps in the 1 Sq 
wave0. This is expected since the proton fractions (and their respective Fermi 
momenta) are rather small, see Fig. 6. 

For neutrons however, the Fermi momenta are larger, and we would expect 
relativistic effects to be important. At Fermi momenta which correspond to 
the saturation point of nuclear matter, kf< = 1.36 fm _1 , the lowest relativistic 
correction to the kinetic energy per particle is of the order of 2 MeV. At den- 
sities higher than the saturation point, relativistic effects should be even more 
important, as can clearly be seen in the calculations of Ref. [72]. Since we are 
dealing with very small proton fractions, a Fermi momentum of kp = 1.36 
fm _1 , would correspond to a total baryonic density ~ 0.09 fm -3 . Thus, at 
larger densities relativistic effects for neutrons should be important. This is 
also reflected in Fig. 24 for the pairing gap in the 3 P2 channel. The relativis- 
tic 3 P2 gap is less than half the corresponding non-relativistic one, and the 
density region is also much smaller. This is mainly due to the inclusion of 
relativistic single-particle energies in the energy denominator of Eq. (68) and 
the normalization factors for the Dirac spinors in the NN interaction. As an 
example, at a neutron Fermi momentum kp = 1.5 fm _1 , the gap has a value of 
0.17 MeV when one uses free single-particle energies and a bare NN potential. 
Including the normalization factors in the NN interaction, but employing free 
single-particle energies reduces the gap to 0.08 MeV. If we employ only DBHF 
single-particle energies and the bare NN interaction, the gap drops from 0.17 
MeV to 0.04 MeV. Thus, the largest effect stems from the change in the single- 
particle energies, although, the combined action of both mechanisms reduce 
the gap from 0.17 MeV to 0.015 at kp = 1.5 fm _1 . The NN interaction in 
the 3 P2 channel depends also strongly on the spin-orbit force, see e.g. Fig. 3.3 
in Ref. [39], and relativistic effects tend to make the NN spin-orbit interac- 
tion from the cu-meson in P- waves more repulsive [39]. This leads to a less 
attractive NN interaction in the 3 P2 channel and a smaller pairing gap. 

The present results can be summarized as follows. 

• The 1 Sq proton gap in /3-stable matter is < 1 MeV, and if polarization 



Even smaller differences are obtained for neutrons in the 1 Sr\ channel. 
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Fig. 23. Proton pairing in /3-stable matter for the 1 Sq partial wave. The 
non-relativistic results are taken from Ref. [108]. 

effects were taken into account [97], it could be further reduced by a factor 
2-3. 

• The 3 P2 gap is also small, of the order of ~ 0.1 MeV in /3-stable matter. If 
relativistic effects are taken into account, it is almost vanishing. However, 
there is quite some uncertainty with the value for this pairing gap for den- 
sities above ~ 0.3 fm -3 due to the fact that the NN interactions are not 
fitted for the corresponding lab energies. 

• Higher partial waves give essentially vanishing pairing gaps in /3-stable mat- 
ter. 

• We have omitted a discussion of hyperon pairing, due to the uncertainties in 
the determination of the hyperon-hyperon interaction. We refer the reader 
here to Ref. [110] for a discussion of these gaps. 

Consequences for cooling histories will be discussed in section 5. 



2.8 Quark matter 



When nuclear matter is compressed to densities so high that the nucleon cores 
substantially overlap, one expects the nucleons to merge and undergo a phase 
transition to chiral symmetric and/or deconfined quark matter. Rephrased in 
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Fig. 24. Neutron pairing in /3-stable matter for the 3 i-2 partial wave. The 
non-relativistic results are taken from Ref. [108]. 

terms of the relevant field excitations, we expect a transition from hadronic 
to quark degrees of freedom at high densities. Knowledge of the EOS of both 
hadronic and quark matter is necessary to estimate the possible effects of this 
transition in neutron stars. 

Recent advances in the QCD phase diagram include improved lattice QCD cal- 
culations, random matrix models [111], and models addressing the possibility 
of color superconductivity at fintie density [112]. Lattice QCD can only treat 
the case of zero baryon chemical potential and is therefore not useful for neu- 
tron stars. Lattice calculations suggest that QCD has a first order transition 
at finite temperature and zero chemical potential, provided that the strange 
quark is sufficiently light [113]. The transition weakens and might change to 
second order for large strange quark masses in the limit of QCD with two 
massless flavors. 

Recent work using chiral random matrix models (chRMM) [111] suggests an 
effective thermodynamic potential of the form 



n(cf>;tJ,,T)/N f = ft o (/i,T) + 2 --ln{[0 2 -G U +*T) 5 



! -(/.-iT) 2 ]}.(69) 



Here Qq(/i,T) is unspecified and independent of the chiral mean field 0. The 
scale of dimensional quantities cannot be determined within the chRMM, and 



55 



must be estimated by relating the variance of the Gaussian random matrix 
ensemble to the vacuum expectation value of the field via the Casher-Banks 
relation [111]. The value of at the minima of f2(0;/i, T) is related to the 
quark condensate, (ipip), which is the order parameter for chiral symmetry 
breaking. Minimization of Eq. (69) leads to a fifth order polynomial equation 
for which is identical in form to the results of Landau-Ginzberg theory 
using a 6 potential. One solution to this equation corresponds to the restored 
symmetry phase with = 0. This model predicts a second order transition 
for n — at a temperature, T c , which is generally agreed to be in the range 
140-170 MeV. For T = 0, a first order transition occurs at some Hq. Since the 
phases in which chiral symmetry is broken and restored must be separated 
in the (//, T) plane by an unbroken line of phase transitions, this implies the 
existence of a tricritical point in the theory of massless quarks. 

There has recently been speculation regarding color superconductivity at medium 
densities resulting from non-perturbative attraction between quarks. At finite 
chemical potential, this invariably leads to the possibility of a diquark con- 
densate which breaks global color invariance [112]. The associated color gap is 
~ 100 MeV and may become the thermodynamically favorable phase at high 
baryon densities. When the strange quark is taken into account many different 
phases may exist [114] and such effects require more analysis. 



2.8.1 Bag Models 

Since we do not have a fully reliable theory for the quark matter phase, we 
will for simplicity employ the simple Bag model in our actual studies of the 
mixed phase and neutron start properties. In the bag model the quarks in 
the hadrons are assumed to be confined to a finite region of space, the so- 
called 'bag', by a vacuum pressure B. The pressure from the quarks inside 
the bag is provided by the Fermi pressure and interactions computed to order 
a s = g 2 /4:7T where g is the QCD coupling constant. The pressure for quarks of 
flavor /, with / = u, d or s is [115] 



1 



4tt 2 

7T 3 



Hfkfdlj — 2.5mj) + l.brrifln 




rrifln 



' Vf + kj 

, m f , 



'Hf + k 



f 



771/ 



K f 



(70) 



The Fermi momentum is kf — y^t 

bag constant B simulating confinement is 
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m 2 f) 1 / 2 . The total pressure, including the 



P = P P 



EPf 

f 



B. 



(71) 



56 



The electron pressure is 

p. = S- (72) 

A Fermi gas of quarks of flavor i has density rii = kp i /n 2 , due to the three color 
states. There is no one-gluon exchange interaction energy between quarks of 
different flavor, while that between quarks of flavour i is given by (2a/37r)Ei 
per quark i [116]. Here Ei is the average kinetic energy per quark, and a is the 
strong interaction coupling constant, assumed to have a value of 0.5. The u and 
d quarks are taken to be massless, and s quarks to have a mass of 150 MeV. 
Typical quark chemical potentials fiq^m^/3 are generally much larger. The 
value of the bag constant B is poorly known, and we present results using 
three representative values, £ 1/4 = 100 MeVfm~ 3 [117], fi 1/4 = 150 MeVfm~ 3 
and B l ' A = 200 MeVfm" 3 [118]. 

Another possible model which has been applied to neutron star studies, and 
which differs from the Bag-model is a massive quark model, the so-called 
Color-Dielectric model (CDM) [119-121]. The CDM is a confinement model 
which has been used with success to study properties of single nucleons, such 
as structure functions [122] and form factors [123], or to describe the inter- 
action potential between two nucleons [124], or to investigate quark matter 
[121,125]. In particular, it is possible, using the same set of parameters, both 
to describe the single nucleon properties and to obtain meaningful results for 
the deconfinement phase transition [121]. The latter happens at a density of 
the order of 2-3 times no when symmetric nuclear matter is considered, and at 
even smaller densities for matter in /9-equilibrium, as discussed below in this 
work. Another important feature is that effective quark masses in the CDM 
are always larger than a value of the order of 100 MeV, hence chiral symmetry 
is broken and the Goldstone bosons are relevant degrees of freedom. This is 
to be contrasted with models like the MIT bag, where quarks have masses 
of a few MeV. We therefore expect the CDM to be relevant for computing 
the cooling rate of neutron stars via the Urea mechanism, as suggested by 
Iwamoto [126]. 

We will however stick to the Bag-model in our discussion of properties of 
neutron stars. 



THERMODYNAMICS OF MULTI-COMPONENT PHASE TRAN- 
SITIONS 



Numerous phase transitions may occur in neutron stars, e.g. the nuclear liquid- 
gas transition in the inner crust and in the interior quark matter and/or con- 
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densates of kaons, pions, hyperons, etc. may be present. Even if the transition 
between two such phases is first order microscopically, the transition macro- 
scopically (in bulk) may be second order when there are several components 
present. We shall in this section briefly describe this curious phenomenon, the 
thermodynamics of multi-component systems and the corresponding mixed 
phases. It will be employed for neutron stars in the subsequent section. 



N H H Q 







Fig. 25. Sketch of a neutron, hybrid, mixed and quark star. 



3.1 Maxwell construction for one component systems 



In the usual picture the transition between two phases occurs at a unique 
pressure, temperature and chemical potential. Consequently the density is 
expected to jump disco ntinuously at the boundary between the two phases. 
This is not only true for systems of one component as in the everyday example 
of water freezing or evaporating. It is also the case for some two component 
systems as, e.g. electrically neutral nuclear matter in /9-equilibrium. Electric 
neutrality requires that the proton and electron densities are the same in bulk 

n p = n e . (73) 

/9-equilibrium requires that the chemical potentials of neutrons and protons 
only differ by that of the electrons 

/i n = /i p + /if M , (74) 

in nuclear matter. These two conditions restrict two of the three components 
leaving only one independent variable, the baryon density. 
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Likewise, in quark matter charge neutrality implies that 
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/3-equilibrium requires analogously 

LL d = Li s = Li u + He M ■ 



(75) 



(76) 



Over the past two decades many authors have considered the properties of 
neutron stars with a core of quark matter [9] . In such "hybrid" stars (see Fig. 
25) it is assumed that each of the two phases are electrically neutral separately 
as in Eqs. (73) and (75) and in /^-equilibrium separately as in Eqs. (74) and 
(76). Gibb's conditions Pnm = Pqm and fi^ M = [i® M (the temperatures are 
vanishing in both phases) then determine a unique density at which the two 
bulk neutral phases coexist. This is the standard Maxwell construction and is 
seen in Fig. 3.1 as the double-tangent. In a gravitational field the denser phase 
(QM) will sink to the center whereas the lighter phase (NM) will float on top 
as a mantle as icebergs in the sea. At the phase transition there is a sharp 
density discontinuity and generally jj,^ m ^ ii® M so that the electron densities 
n e = hI/3ti 2 are different in the two phases. This assumes that the sizes of 
QM structures are larger than electron screening lengths which, as discussed 
in [10], is not the case. 

For small values of the bag constants the phase transition occurs at densities 
lower than no and the whole neutron star is a quark star except possibly for 
a hadronic crust [127]. 
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Fig. 26. Pressure as function of density. Maxwell construction and mixed phase 
for bag parameter B 1 ' = 200 MeVfm" 3 . The arrow indicates where the Maxwell 
construction starts. 
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3.2 Two component systems in a mixed phase 



For several coexisting components as, e.g. dissolved chemicals, there may ex- 
ist a mixed phase where the various chemical potentials of the solvents vary 
continuously - as does the pressure and densities. A similar phenomenon was 
predicted for the nuclear liquid-gas phase transitions in the inner crust of neu- 
tron stars confirmed by recent detailed numerical calculations [33,34]. In the 
inner crust the nuclei are surrounded by a neutron gas with an interpene- 
trating constant background of electrons, i.e. the nuclear matter and neutron 
gas form a mixed phase. Going a few hundred meters down in the neutron 
star crust the density of nuclei and thus also the average density increase 
continuously. Thus there is no sharp density discontinuity in bulk, i.e. over 
macroscopical distances of more than hundreds of Fermi's and up to several 
meters (see Fig. 25). However, on microscopical distances of a few Fermi's the 
density varies rapidly. Another example first considered by Glendenning [5] is 
that of a mixed phase of nuclear and quark matter. 

Contrary to the Maxwell construction described in the previous subsection, 
where the condition of charge neutrality applies to both phases, Eqs. (73) and 
(75), it is relaxed to overall charge neutrality only. Thus two conditions are 
relaxed to one - allowing for one new variable quantity, which is usually taken 
as the filling fraction f of one of the phases in their coexisting mixture. For 
the nuclear and quark matter mixed phase the filling fraction is defined as the 
fraction of the volume which is in the quark phase 

Z-F^T^' (77) 

and so (1 — /) is the filling fraction of nuclear matter. The overall charge 
neutrality requirement is 

2 11 

fn p + (1 - /) {-n u - -n d - ~n s ) = n e . (78) 

A number of requirements must be met in order to form such a mixed phase 
as addressed in [10]. These will be discussed in the following section for the 
mixed phase of nuclear and quark matter. 



4 STRUCTURE OF NEUTRON STARS 



The structure of the neutron star is seriously affected by phase transitions, 
the order of the phase transition and whether mixed phases can occur over 
a significant part of the star. The important questions to be addressed are 
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numerous. When is it legitimate to regard the electron density as uniform, 
what is the spatial structure of the new phase, and is it energetically favorable? 
In order to answer these questions one must first investigate screening lengths 
of the various charged particles and compare to typical size scales of structures. 
For that matter Coulomb and surface energies must be calculated. 

We will mainly discuss the mixed phase of nuclear and quark matter in cores 
of neutron stars. The mixed phase of nuclei and a neutron gas is in many ways 
similar and has been calculated in detail in Refs. [33,34,129] 



4-1 Screening lengths 



As described in [10] the mixed phase of quark and nuclear matter may be 
regarded as droplets of quark matter immersed in nuclear matter at lower 
densities, usually referred to as the droplet phase, even though at higher den- 
sities its structure is more complicated. If droplet sizes and separations are 
small compared with Debye screening lengths, the electron density will be 
uniform to a good approximation. 

The Debye screening length, A^ is given by 

where Hi, fii, and Qi are the number density, chemical potential, and charge 
of particle species %. Considering only electrons gives a screening length 

Ag> = ±J— , (80) 

K F,e 

where a ~ 1/137 and the Fermi momentum kF,e — He since the electrons are 
always relativistic at these densities. For /i e £150 MeV we thus obtain A^ £13 
fm. The screening length for protons alone Ajy , is given by (irvF, p /c4:a(l + 
Fq)) 1 ' 2 /kF,p, where Fq is the Landau parameter which gives the energy for 
proton density variations. At the saturation density for symmetric nuclear 
matter, Fq ~ 0, whereas at higher densities Fq ~ 1 [128]. Since \x v ~ m, 

(v) 
the nucleon mass, we find A^, £10 fm, somewhat shorter than the electron 

screening length, and therefore in the nuclear matter phase, protons are the 

particles most effective at screening. The screening length for quarks is \£ ~ 

7/kF, q where q—u, d, and s refer to up, down and strange quarks. It depends 

only slightly on whether or not s-quarks are present, so for fi q ~ m/3 we find 

A§ } ~ 5 fm. 
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In a composite system, such as the one we consider, screening cannot be de- 
scribed using a single screening length, but it is clear from our estimates, that 
if the characteristic spatial scales of structures are less than about 10 fm for 
the nuclear phase, and less than about 5 fm for the quark phase, screening ef- 
fects will be unimportant, and the electron density will be essentially uniform. 
In the opposite case, when screening lengths are short compared with spatial 
scales, the total charge densities in bulk nuclear matter and quark matter will 
both vanish. 



4-2 Surface and Coulomb energies of the mixed phase 



When screening lengths are much larger than the spatial scale of structures. 
This condition implies that the electron density is uniform everywhere, and 
all other particle densities are uniform within a given phase. The problem 
is essentially identical to that of matter at subnuclear densities [129], and 
the structure is determined by competition between Coulomb and interface 
energies. When quark matter occupies a small fraction, /, of the total volume, 
it will form spherical droplets immersed in nuclear matter. For higher filling 
fractions, the quark matter will adopt shapes more like rods ( "spaghetti" ) and 
plates ("lasagna"), rather than spheres. For / > 0.5, the structures expected 
are the same as for a filling factor 1 — /, but with the roles of nuclear matter 
and quark matter reversed. Thus one expects for increasing / that there will 
be regions with nuclear matter in rod-like structures, and roughly spherical 
droplets. 

To estimate characteristic dimensions, some special cases can be considered. 
The intricate structures in the general case will be discussed in Sec. 4.5. When 
/ is small or close to unity, the minority phase will form spherical droplets. 
The surface energy per droplet is given by 

S s = o-AnR 2 , (81) 

where a is the surface tension, and the Coulomb energy is 

3 Z 2 e 2 16tt 2 2 5 

£c = t— B - = ~TF~ \PQM - Pnm) R ■ (82) 

Here Z is the excess charge of the droplet compared with the surrounding 
medium, Ze = (pqm — Pnm)Vd where Vd = (47r/3)i? 3 is the droplet volume 
and pqm and p NM are the total charge densities in bulk quark and nuclear 
matter, respectively. Minimizing the energy density with respect to R one 
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obtain the usual result that Ss = 2£c and find a droplet radius 



R 



15 a \ 1/3 

8tt(p qm - Pnm) 2 ) 



^5.0fm^V /3 f PQM ' PJVM V 2/3 . (83) 

Wo/ V Po / 

In the second formula we have introduced the quantities po — e 0.4 fm -3 
and o"o = 50 MeV-fm -2 which, as we shall argue below, are typical scales 
for the quantities. (A droplet of symmetric nuclear matter in vacuum has 
a surface tension a = 1 MeV-fm -2 for which (83) gives R ~ 4 fm, which 
agrees with the fact that nuclei like 56 Fe are the most stable form of matter 
for roughly symmetric nuclear matter at low density.) The form of Eq. (83) 
reflects the fact that on dimensional grounds, the characteristic length scale 
is (ct/(pq M — pnm) 2 ) 1 ^ 3 times a function of /. The total Coulomb and surface 
energy per unit volume is given for small / by 



£s+c = f 9 (—ct 2 (pqm ~ Pnm) 2 ) 



1/3 

2/3 



~ 44MeV fm" 3 / ^ Pqm ~ PjvM I . (84) 

The result for / close to unity is given by replacing / by 1 — /. In the case when 
the volumes of quark and nuclear matter are equal, / = 1/2, the structure 
can be approximated as alternating layers, of quark and nuclear matter and 
was considered in [10]. 

To estimate length scales and energy densities, one needs the surface tension of 
quark matter and the charge densities in the two phases. A rough estimate of 
the surface tension is the bag constant, B, times a typical hadronic length scale 
~1 fm. Estimates of the bag constant range from 50 to 450 MeV-fm -3 [130]. 
The kinetic contribution to the surface tension at zero temperature has been 
calculated in the bag model in Ref. [131]. Only massive quarks contribute 
because relativistic particles, unlike non-relativistic ones, are not excluded 
near the surface due to the boundary conditions. The kinetic contribution to 
a from a quark species depends strongly on its mass and chemical potential. 
For m s <C p s it behaves as (3/47i 2 )p 2 s m s , and it vanishes as m s approaches 
p s . If we adopt for the strange quark mass the value m s ~ 150 MeV, and 
for the quark chemical potentials one third of the baryon chemical potential, 
which generally is slightly larger than the nucleon mass, p s ~ Pb/3 £ to/ 3, 
one obtains from Ref. [131] o ~ 10 MeV-fm -2 , which is close to the maximum 
value it can attain for any choice of m s . We conclude that the surface tension 
for quark matter is poorly known, but lies most probably in the range 10 — 100 
MeV-fm -2 . (Lattice gauge theory estimates of a at high temperatures and 
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zero quark chemical potentials lie in the range a ~ 0.14 — 0.28T C 3 ~ 10-60 
MeV-fm~ 2 [132] for T c ~ 150 — 200 MeV, comparable to our estimates for cold 
quark matter, but it is unclear to what extent this agreement is accidental.) 

To estimate charge densities, we consider quark matter immersed in a uniform 
background of electrons, /^-equilibrium insures that \Xd — p, s — p, u + p, e , and 
therefore in the absence of quark-quark interactions, one finds the total electric 
charge density in the quark matter phase is given for \x e <C p u ~ pd = Pq and 
m s < p q by 



Pqm = -(2n u -n d -n s - 3n e ) ~ — ( -m\p q - 2p e p 2 ) . (85) 



e e /l 

— {2n u — rid — n s — 3n e ) ~ — \-..„ sr , q -r~ ct ~ q 

Assuming m s ~ 150 MeV and ii q ~ m/3 the second term dominates except 
for small ii e and so the droplet is negatively charged and for ii e ~ 170 MeV 
the density is about — 0.4e fm~ 3 , the characteristic scale of densities adopted 
in making estimates above. 

Due to the high quark density, p^M is small compared with pq M in Eq. (83) 
when quark matter occupies a small fraction of the volume. The electron 
chemical potential in neutron stars depends strongly on the model for the 
nuclear equation of state, but generally one finds /x e ^170 MeV. Consequently, 
for a ~ 10 MeV-fm -2 we find from Eq. (83) a radius of R^ 3.1 fm, whereas 
a ~ 100 MeV-fm -2 gives R^, 6.6 fm. For / close to unity one finds nuclear 
bubble radii which are comparable with those for quark droplets, and for 
the layer-like structures expected for / ~ 0.5, half the layer thickness is of 
comparable size. Estimates of characteristic scales for rod-like structures give 
similar values. 

Detailed calculations show that the effects of nonuniformity of the charge 
distribution affect estimates of Coulomb energies significantly if the charac- 
teristic lengths, R and a, exceed the Debye screening length. The estimates 
of screening lengths made above show that screening will be not be dominant 
for surface tensions below about 100 MeV, if the charge density difference is 
Po, but for higher values the simple picture of coexisting uniform bulk phases 
would become invalid, and the droplet phase would resemble increasingly two 
electrically neutral phases in equilibrium. 

For the smallest droplets of size R ~ 5fm the charge Z = (4n/3)R 3 pQM and 
baryon number A = (47r/3)R 3 n are typically a few hundreds. 
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Jf.,3 Is the mixed phase energetically favored? 



If the bulk energy gained by going to the mixed phase is larger than the 
costs of the associated Coulomb and surface energies, then the mixed phase 
is energetically favored. Before going on to estimate these crucial energies 
we point out the basic physical reason why the bulk energy is lower in the 
mixed phase. As mentioned in connection with Eq. (85) the QM droplets 
are negatively charged. By immersing QM in the positively charged NM we 
can either remove some of the electrons from the top of the Fermi levels with 
energy // e , or we can increase the proton fraction in NM by which the symmetry 
energy is lowered. In equilibrium a combination of both will occur and in both 
cases bulk energy is saved and a lower energy density is achieved as seen in 
Fig. 3.1. 

To calculate the bulk energy we adopt a simple form for the energy density of 
nuclear matter consisting of a compressional term, a symmetry term, and an 
electron energy density as discussed in Eq. (49) 



£nm = nE = n[m + E comp (n) + S(n)x 2 ] + e 



comp\ 

4 




= n[m + £„_(_ - 2 - S)(l + 5-)- L + 5 (-)V] + -^. (86) 

The electron chemical potential is never much above the muon mass and there- 
fore muons may be ignored. For quark matter we assume the bag model equa- 
tion of state 

e Q M = {l~— )| £ SI + 5 + 72^> ( 87 ) 

with the QCD fine structure constant a s ~ 0.4 and bag constant B 1 ^ 4 ~ 120 
MeV-fm -3 . We have taken all quark masses to be zero. In the absence of 
surface and Coulomb effects the equilibrium conditions for the droplet phase 
are that the quark and nuclear matter should have equal pressures, and that 
it should cost no energy to convert a neutron or a proton in nuclear matter 
into quarks in quark matter. The last condition amounts to fi n = 2/i^ + /i M and 
Up = Hd + 2/i n - The electron density is the same in quark and nuclear matter, 
and we assume that matter is electrically neutral and in /9-equilibrium, that 
is ji n = jj,p + jj, e and fid = fi u + /x e . The chemical potentials are related to 
the Fermi momenta by \x q = PF,q(^ — 2a s /7r) -1 / 3 . Electrons contribute little 
to pressures, but they play an important role through the /9-equilibrium and 
charge neutrality conditions. 

Figure 27(a) shows the density dependence of the energy density of the droplet 
phase calculated neglecting surface and Coulomb energies (a = 0) for a simple 
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Fig. 27. (a) The full line gives the energy density of the droplet phase without 
surface and Coulomb energies (a = 0). Also shown are the energy densities of 
electrically neutral bulk nuclear matter, quark matter in /3-equilibrium, and the 
double tangent construction (dashed line) corresponding to the coexistence of bulk 
electrically neutral phases, (b) Energy densities of the droplet phase relative to its 
value for a= for a = 10, 50, and 90 MeV-fm~ 2 . When the energy density of the 
droplet phase falls within the hatched area it is energetically favored. From [10]. 

quadratic EoS for nuclear matter [10]. The energy of uniform, electrically 
neutral, bulk nuclear matter in /9-equilibrium is also shown, together with the 
corresponding result for quark matter. The double-tangent construction gives 
the energy density for densities at which the two bulk neutral phases coexist. 
This corresponds to the standard treatment of the phase transition between 
nuclear matter and quark matter, in which the pressure remains constant 
throughout the transition, and consequently neutron stars have a core of quark 
matter and a mantle of nuclear matter, with a sharp density discontinuity 
at the phase transition. As one sees, if surface and Coulomb effects may be 
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ignored, the transition from nuclear matter to the droplet phase occurs at a 
lower density than the transition to two bulk neutral phases, a feature also 
apparent in Ref. [5]. In addition, droplets of nuclear matter survive up to 
densities above those at which bulk neutral phases can coexist. We also observe 
that bulk contributions to the energy density of the droplet phase are always 
lower than those for coexisting bulk neutral phases. While detailed properties 
of the droplet phase depend strongly on the bulk energies, the qualitative 
picture we find persists over a wide range of possible bulk matter properties. 

We now estimate surface and Coulomb energies. When quark matter occupies 
a small fraction of space, /, one can show that the difference in energy between 
the droplet phase and bulk neutral nuclear matter varies as f 2 . In contrast 
to this, the contributions to the energy density from surface and Coulomb 
energies are linear in /, see Eq. (84). Similar results apply for / close to unity. 
This shows that the transitions to the droplet phase must occur via a first- 
order transition. However, if the surface and Coulomb energies are sufficiently 
large, the droplet phase may never be favorable. The energy-density difference 
between the droplet phase, neglecting surface and Coulomb effects, and two 
coexisting neutral phases is at most 10 MeV-fm -3 , as may be seen from Fig. 
27. This is very small compared with characteristic energy densities which 
are of order 1000 MeV-fm -3 . In Fig. 27(b) we show the energy density of the 
droplet phase for various values of the surface tension, relative to the value for 
a = 0. In these calculations the geometry of the droplets was characterized by 
a continuous dimensionality, d, as described in Ref. [129], with d = 3, 2 and 1 
corresponding to spheres, rods and plates, respectively. For the droplet phase 
to be favorable, its energy density must lie below those of nuclear matter, 
quark matter, and coexisting electrically neutral phases of nuclear and quark 
matter. That is the droplet phase will be favored if its energy lies within the 
hatched region in Fig. 27(a+b). We see that whether or not the droplet phase 
is energetically favorable depends crucially on properties of quark matter and 
nuclear matter. For our model the droplet phase is energetically favorable at 
some densities provided <r^,70 MeV-fm -2 . However, given the large uncertain- 
ties in estimates of bulk and surface properties one cannot at present claim 
that the droplet phase is definitely favored energetically. 



4-4 Melting temperatures 



Should the quark-droplet phase exist in neutron stars, it could have impor- 
tant observational consequences. First, as Glendenning showed, the pressure 
difference across the droplet phase can be large, of order 250 MeV-fm~ 3 . This 
is also seen from Fig. 27(a), since the pressure is the negative intercept of the 
tangent to the curve. Consequently, a large portion of a neutron star could 
consist of matter in the droplet phase. Secondly, phases with isolated droplets 
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would be expected to be solid. The melting temperature is ~ Z 2 e 2 f 1 ^ 3 /(170R) 
[133] typically some hundreds of MeV, while spaghetti- and lasagna-like struc- 
tures would exhibit anisotropic elastic properties, being rigid to some shear 
strains but not others in much the same way as liquid crystals. This could be 
important for quake phenomena, which have been invoked to explain observa- 
tions in a number of different contexts. Third, neutrino generation, and hence 
cooling of neutron stars could be influenced. This could come about because 
nuclear matter in the droplet phase has a higher proton concentration than 
bulk, neutral nuclear matter, and this could make it easier to attain the thresh- 
old condition for the nucleon direct Urea process [134]. Another is that the 
presence of the spatial structure of the droplet phase might allow processes to 
occur which would be forbidden in a translationally invariant system. Finally, 
one should bear in mind the possibility that even if the droplet phase were 
favored energetically, it would not be realized in practice if the time required 
to nucleate is too long. 



4-5 Funny phases 



Surface and Coulomb energies determine the topology and length scales of 
the structures. Denoting the dimensionality of the structures by d (d — 3 for 
droplets and bubbles, d = 2 for rods and d — 1 for plates) the surface and 
Coulomb energies are generally [26,33] 
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where a is the surface tension, R the size of the structure, and pqm and p^M 
are the total charge densities in bulk QM and NM, respectively. For droplets 
(/ ~ 0) or bubbles (/ ~ 1) d = 3 and the Coulomb energies reduce to the usual 
term Eq = (3/5)Z 2 e 2 /R where Z is the excess charge of the droplet compared 
with the surrounding medium, Ze = (4tt/3)(pqm — pnm)R 3 - Minimizing the 
energy density with respect to R we obtain the usual result that £s = 2£c- 
[Q Minimizing with respect to the continuous dimensionality as well thus 
determines both R and d 

For the Walecka model the droplet phase is energetically favorable at some den- 
sities provided <t^20 MeV/fm 2 [135]. For comparison, using a quadratic EoS 
for NM[10] one finds in stead the more favorable condition a ^70 MeV/fm 2 . 
Given the large uncertainties in estimates of bulk and surface properties one 



The condition for fission instability is contrarily: 2£$ < £c 
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cannot at present claim that the droplet phase is definitely favored energeti- 
cally. 

The mixed phase in the inner crust of neutron stars consists of nuclear matter 
and a neutron gas in /^-equilibrium with a background of electrons such that 
the matter is overall electrically neutral [33,34]. Likewise, quark and nuclear 
matter can have a mixed phase [5] and possible also nuclear matter with and 
without condensate of any negatively charged particles such as K~ [91], n~ , 
E~, etc. The quarks are confined in droplet, rod- and plate-like structures [10] 
(see Fig. 28) analogous to the nuclear matter and neutron gas structures in 
the inner crust of neutron stars [33,34]. Depending on the equation of state, 
normal nuclear matter exists only at moderate densities, p ~ 1 — 2p - With 
increasing density, droplets of quark matter form in nuclear matter and may 
merge into rod- and later plate-like structures. At even higher densities the 
structures invert forming plates, rods and droplets of nuclear matter in quark 
matter. Finally, pure quark matter is formed at very high densities unless the 
star already has exceeded its maximum mass. 

A necessary condition for forming these structures and the mixed phase is that 
the additional surface and Coulomb energies of these structures are sufficiently 
small. Excluding them makes the mixed phase energetically favored [5]. That 
is also the case when surface energies are small (see [10] for a quantitative 
condition). If they are too large the neutron star will have a core of pure quark 
matter with a mantle of nuclear matter surrounding and the two phases are 
coexisting by an ordinary first order phase transition. 

The quark and nuclear matter mixed phase has continuous pressures and den- 
sities [5] when surface and Coulomb energies are excluded. There are at most 
two second order phase transitions. Namely, at a lower density, where quark 
matter first appears in nuclear matter, and at a very high density, where all 
nucleons are finally dissolved into quark matter, if the star is gravitation- 
ally stable at such high central densities. However, due to the finite Coulomb 
and surface energies associated with forming these structures, the transitions 
change from second to first order at each topological change in structure [10]. 
If the surface and Coulomb energies are very small the transitions will be only 
weakly first order but there may be several of them. 



4-6 Summary of neutron star structures 



To summarize, we have shown that whether or not a droplet phase consisting 
of quark matter and nuclear matter can exist in neutron stars depends not 
only on bulk properties, but also on the surface tension. In order to make 
better estimates it is important to improve our understanding of the transition 
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between bulk nuclear matter and bulk quark matter. For the droplet phase 
to be possible, this must be first-order. If the transition is indeed first-order, 
better estimates of the surface tension are needed to determine whether the 
droplet phase is favored energetically. 

In these analyses several restrictions were made: the interfaces were sharp, the 
charge densities constant in both NM and QM and the background electron 
density was also assumed constant. Relaxing these restrictions generally allow 
the system to minimize its energy further. Constant charge densities may 
be a good approximation when screening lengths are much larger than spatial 
length scales of structures but since they are only slightly larger [10] the system 
may save significant energy by rearranging the charges. 
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Fig. 28. Nuclear and quark matter structures in a ~ 1.4M neutron star. Typical 
sizes of structures are ~ 10 -14 m but have been scaled up to be seen. 



5 OBSERVATIONAL CONSEQUENCES FOR NEUTRON STARS 



In this section we first briefly review the observational status of neutron star 
masses from binary pulsars and X-ray binaries. Subsequently, we revisit the 
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equations for calculating the structure of rotating neutron stars and calculate 
masses, radii, moments of inertia, etc. for rotating neutron stars with the 
equations of state described in the previous sections with and without phase 
transitions. Glitches are then discussed in subsections 5.6 and 5.7 while a 
discussion on neutron star cooling is given in subsection 5.8. In the last two 
subsections we discuss supernovae and Gamma-Ray-Bursters. 



5. 1 Masses from radio pulsars, X-ray binaries and QPO 's 



The measurements of masses and radii of neutron stars (as well as detailed 
study of their cooling histories and rotational instabilities) may provide a 
unique window on the behavior of matter at densities well above that found 
in atomic nuclei. The most precisely measured physical parameter of any rotat- 
ing neutron star, or pulsar, is its spin frequency. The frequencies of the fastest 
observed pulsars (PSR B1937+21 at 641.9 Hz and B1957+20 at 622.1 Hz) 
have already been used to set constraints on the nuclear equation of state at 
high densities under the assumption that these pulsars are near their maxi- 
mum (breakup) spin frequency. However, the fastest observed spin frequencies 
may be limited by complex accretion physics rather than fundamental nuclear 
and gravitational physics. A quantity more directly useful for comparison with 
physical theories is the neutron star mass. In Fig. 5.1 we show the latest com- 
pilation of Thorsett and Chakrabarty [136] of neutron star masses in binary 
radio pulsar systems. As can be seen from this figure, most of observed binary 
pulsars exhibit masses around 1.4M . One exception may be PSR J1012+5307 
with mass M = (2.1 ±O.4)M [137]. However, its mass is less well determined 
as the system is non-relativistic. The masses are determined from pulse delays 
of the millisecond pulsar as well as radial- velocity curves and spectral lines of 
the white dwarf companion 

The recent discovery of high-frequency brightness oscillations from 11 neutron 
stars in low-mass X-ray binaries may provide us with a new promising method 
for determining masses and radii of neutron stars. These quasi-periodic oscilla- 
tions (QPO) are observed in both the persistent X-ray emission and in bursts. 
According to the most successful model of Miller et al. [138] the QPO's are 
most likely the orbital frequencies of accreting gas in Keplerian orbits around 
neutron stars. The orbital frequency of the gas at distance r from the neutron 
star is 

, Q po = ^(-) ■ (90) 

As the QPO last many periods, the gas has to be in a stable orbit. The 
innermost stable orbit R ms for a slowly rotating neutron star is related to its 
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Fig. 29. Neutron star masses from observations of radio pulsar systems. Error bars 
indicate central 68% confidence limits, except upper limits are one-sided 95% con- 
fidence limits. The vertical lines are drawn at a mass 1.35 ± O.O4M0. Taken from 
Ref. [136]. 



mass as 



M — —R m s- 



(91) 



Thus we obtain the limits on the non-rotating neutron star masses and radii: 



kHz 



M < 2.2M G 


^QPO 




kHz 


R < 20km- 





^QPO 



(92) 
(93) 



For example, the 1220 Hz QPO observed in the atoll source 4U 1636-536 limits 
the neutron star mass to M < 1.8M and its radius to R < 16km. It was 
predicted by Miller et al. [138] that as the accretion flux increases towards 
the innermost stable orbit the QPO frequency should stop increasing. This 
was subsequently observed by Zhang et al. [3] and Kaaret et al. [139] for 4U 
1820-30 which has vqpo = 0.8 — 0.9 kHz. The resulting neutron star mass is 
M ~ 2.3M when rotation is included. The deduced mass is consistent with 
the hypothesis that these neutron stars were born with M ~ 1.4M Q and have 
been accreting matter at a fraction of the Eddington limit for 10 8 yr. The 
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QPO's provide an important tool for determining neutron star masses or at 
least restricting them and for limiting radii. Due to accretion they are expected 
to be heavier than X-ray binaries and therefore potentially more interesting 
for restricting maximum masses and the EoS. 

A new determination of the mass of the accretion powered X-ray pulsar Vela 
X-l gives the mass (95% confidence limit) 

M = 1.93+°;^M . (94) 

The mass could be determined from the velocities of the binaries from pulse 
timing of the neutron star and doppler shifts in the spectral lines of its com- 
panion as well as the inclination from orbital flux variations. 

Recently, also the mass of the low mass X-ray binary Cygnus X-2 has been 
estimated from U-B-V light curves [140] 

M = (1.78 ± 0.23) M Q . (95) 

The existence of such large neutron star masses will require a rather stiff EoS 
for nuclear matter and restricts the softening due to phase transitions severely, 
as will be discussed in the following subsection. 



5. 2 TO V and Hartle 's equations 



The theoretical description of a neutron star is governed by conditions imposed 
by general relativity. General relativity has to be taken into account for the 
determination of the gross properties of a star with approximately one solar 
mass Mq and a radius R of approximately 10 km, since relativistic effects are 
of the order [46] 

M 

— -0.1-0.2. (96) 

R 

The starting point for such studies is how to determine Einstein's curvature 
tensor G^ for a massive star (it^„, g^,,, and R denote the Ricci tensor, metric 
tensor, and Ricci scalar, respectively). 

GV = R^ - - 9llu R = 8tt TV (e, P(e)) . (97) 

A necessary ingredient for solving this equation is the energy-momentum ten- 
sor density T^,, for which knowledge of the EoS, i.e. pressure P as function of 
the energy density e is necessary. 
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For a spherically symmetric and static star, the metric has the Schwarzschild 
form 

ds 2 = _ e W dt 2 + e 2A(r) rfr 2 + ^2 + ^2^2) ^ (gg) 

where the metric functions are given by: 

e 2A W = (l- 7 (r))- 1 , (99) 

e 2^) =e -2A(r) = (1 _ 7(r)) for r>Rstar) (100 ) 

with 

f 2M(r)/r , r < R s 
7 (r) = (101) 

[ 2M s /r ,r>R s 

Einstein's equations for a static star reduce then to the familiar Tolman- 
Oppenheimer-Volkoff equation (TOV) [46,141,142]: 

^^ = -\ (e(r) + P(r)) (Mir) + 4nr 3 P(r)) e 2A(r) , (102) 

dr r v ' 

where the gravitational mass M(r) contained in a sphere with radius r is 
determined via the energy- density e(r) by: 

r 

M(r) = 4tt f e{r)r 2 dr . (103) 



The metric function 0(r) obeys the differential equation 

d<p _ 1 dP 

dr e{r) + P(r) dr 

with the boundary condition 



(104) 



(f>(r = R s ) = ^]n^-'Y(Rs))- (105) 

For a given EoS, i.e. P(e), one can now solve the TOV equation by integrating 
them for a given central energy density e c from the star's centre to the star's 
radius, defined by P(R S ) = 0. 

More complicated is the case of rotating stars, where due to the rotation 
changes occur in the pressure, energy density, etc. The energy-momentum 
density tensor T^ takes the form {g^u^Uv = — 1) [127,143,144]: 

T^ = T° + AT^, (106) 
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with 



7%, =(e + P) UfM u v + P 9ilu , 
AT M , = (Ae + APKu„ + AP<? 



(107) 
(108) 



P, e, and p are quantities in a local inertial frame comoving with the fluid at 
the instant of measurement. For the rotationally deformed, axially-symmetric 
configurations one assumes a multipole expansion up to second order (P 2 de- 
notes the Legrendre polynomial): 



AP = {e + P){ Po +p 2 P 2 {co S 6)) 



(109) 
(110) 

(111) 



For the rotating and deformed star with the rotational frequency Q one has 
now to deal with a generalized Schwarzschild metric, given by [145,146] 






(112) 



Here, cj(r) denotes the angular velocity of the local inertial frame, which, due 
to the dragging of the local system is proportional to Q. 

The metric functions of Eq. (112) which correspond to stationary rotation and 
axial symmetry with respect to the axis of rotation are expanded up to second 
order as (independent of <f> and t) 



,Mrm = e 20M [! + 2 (/lo ( r> n) + h ^ fi )p 2 (cos0))] , 
r 2 sin 2 6 [1 + 2 (v 2 (r, H) - h 2 (r, ft)) P 2 (cos 0)] 



3 2</>(r,<j!>,n) 



; 2Mr,e,n) = r 2 [! + 2 ( U2 ( r> ft) _ /i 2 ( r , ft)) P 2 (cos( 



3 2A(r,0,Q) 



3 2A(r) 



2 mo(r, ft)G + m2(r, ft)P 2 (cos i 
r 1 — 7(r) 



(113) 
(114) 
(115) 

(116) 



The angular velocity in the local inertial frame is determined by the differential 
equation 



^l r 4 j(r )^M 4r 3^ 



d 

dr 



dcu\ 

dr J 



dr 



u(r) = , r < R s 



(117) 



where u(r) is regular for r = with ^ = 0. j(r) abbreviates 



j[r) =e 



'Vl- 



7(r) . 



'118) 
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Outside the star w(r, Q) is given by: 

2 
uj(r,n) = n -J(ft) , r>R s . (119) 

The total angular momentum is defined by: 



J(n >=f * „ ■ (120) 

\ / r=R s 

From the last two equations one obtains then an angular frequency Q as a 
function of central angular velocity u c = u(r = 0) (starting value for the 
iteration): 

Q(u> c ) = lu(R s ) + J^(^) • (121) 

Due to the linearity of Eq. (117) for oo(r) new values for uj{r) emerge simply 
by rescaling of uj c . The momentum of inertia, defined by I = ^, is given by 

■-^'-^L * + P "-"-^ (122) 



Q 3 J ^1 - 7 (r) ^ 

Relativistic changes from the Newtonian value are caused by the dragging 
of the local systems, i.e. u/Q, the redshift (e _<?i ), and the space-curvature 
( (1 — 7(r))~ ' J. For slowly rotating stars with low masses, one can neglect 
the dragging (^ — ► 1) and rotational deformations, but we would like to 
emphasize that the described treatment is not restricted by low masses and/or 
slow rotations. 

With u(r), one can also solve the coupled mass monopole equations (£ = 0) for 
rrio,po, where the latter represents the monopole pressure perturbation, and 
ho, Refs. [127,143]). The quadrupole distortions Y12 and i>2 [i = 2) determine 
the star's shape (see Refs. [127,143]). We will here just state the equations for 
the monopole functions itlq and po i n order to obtain the corrections to the 
mass due to rotation. We will not deal with quadrupole corrections in this 
work. The equations read 

drriQ A 2 de 1 .,, , 



and 




dr 47r V (£ + P)P0+ 12 Jr ( ~' + - r ^V^7^ ( 123 ) 



dp 1 + 8tt 2 P (g + P) 1 jV fdwV 1 d /rW\ 
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The boundary conditions are itlq — > and po — »■ when r — ► 0. Outside the 
star one has 

m = AM - —J(tt) 2 r>R s . 

with AM being the rotational correction to the gravitational mass. This cor- 
rections is given by 

AM = m + — J(fi) 2 , (125) 

at the surface of the star. Thus, when we solve the monopole equations we 
know also the correction to the gravitational mass. These two equations, to- 
gether with Eqs. (102), (103), (104) and (117) form the starting point for our 
numerical procedure for obtaining the total mass, radius, moment of Inertia 
and rotational mass. Results for various approaches to the EoS are discussed 
in the next subsection. 



5. 3 Neutron star properties from various equations of state 



For a given EoS we obtain the mass and radii of the neutron star by solving 
the equations for a weakly rotating neutron star as given by Hartle [141] and 
discussed in Eqs. (102), (103), (104), (117), (123) and (124). Various results 
with and without rotational corrections are displayed in Figs. 30-38, where we 
show total masses, mass-radius relations and moments of inertia for various 
approximations to the EoS. The following possible properties pertain to the 
various approximations to the EoS. 

• For the EoS parametrization of Akmal et al. [13] with just pn degrees of 
freedom, the EoS with 8 = 0.13 gives the stiff est EoS and thereby the 
largest neutron star mass. For 5 < 0.13 the EoS is superluminal. See also 
the discussion in connection with Figs. 18 and 19. For S = 0.3 or 5 = 0.4 the 
EoS differs from that with S = 0.2 or 5 = 0.13 at densities below n = 0.3 
fm~ 3 . This explains the differences in masses seen at low central densities 
in Fig. 30. The reader should also recall that in section 2, the best fit to the 
results of Ref. [13] was obtained with 5 = 0.2. 

• We have selected three representative values for the Bag-model parameter 
B, namely, 100, 150 and 200 MeVfm" 3 . For B 1 ^ = 100 MeVfm" 3 , the mixed 
phase starts already at 0.22 fm~ 3 and the pure quark phase starts at 1.54 
fm -3 . For B l / A = 150 MeVfm -3 , the mixed phase begins at 0.51 fm -3 and 
the pure quark matter phase begins at 1.89 fm~ 3 . Finally, for B 1 ^ = 200 
MeVfm" 3 , the mixed phase starts at 0.72 MeVfm -3 while the pure quark 
phase starts at 2.11 fm -3 . 

• In case of a Maxwell construction, in order to link the hadronic and the 
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quark matter EoS, we obtain for B 1 ! 4 = 100 MeVfm" 3 that the pure pn 
phase ends at 0.58 fm~ 3 and that the pure quark phase starts at 0.67 fm -3 . 
For B 1 / 4 = 150 MeVfm" 3 , the numbers are 0.92 fm -3 and 1.215 fm~ 3 , 
while the corresponding numbers for B 1 ! 4 = 200 MeVfm~ 3 are 1.04 and 
1.57 fm" 3 . 



As can be seen from Figs. 30, 32, 34, 36 and 38, none of the equations of state 
from either the pure pn phase or with a mixed phase construction with quark 
degrees of freedom, result in stable configurations for densities above ~ 10no, 
implying thereby, see e.g. Figs. 32 and 36, that none of the stars have cores 
with a pure quark phase. The EoS with pn degrees of freedom only results 
in the largest mass ~ 2.2M when the rotational correction of Eq. (125) is 
accounted for, see Fig. 36. With the inclusion of the mixed phase, the total 
mass is reduced since the EoS is softer. 

Several interesting conclusions can be inferred from the results displayed in 
Figs. 30-38. Firstly, to obtain neutron star masses of the order M ~ 2.2M 
as may now have been observed in QPO's, we need the stiff est EoS allowed 
by causality (i.e. 8 ~ 0.13 — 0.2) and to include rotation, see Figs. 30 and 36. 
Furthermore, a phase transition to quark matter below densities of order ~ 5no 
can be excluded, corresponding to restricting the Bag constant to B 1 ^ 4 ^ J 200 
MeVfm -3 . This can be seen in Fig. 32 where we plot star masses as function 
of the central density n c and bag- model parameter B. These results differ 
significantly from those of Akmal et al. and Kalogera and Baym [13,83] due to 
the very different recipes we use to incorporate causality at high densities. In 
Refs. [13,83] the EoS is discontinuously stiffened by taking v s — c at densities 
above a certain value n c which, however, is lower than n s = 5no where their 
nuclear EoS becomes superluminal. This stiffens the nuclear EoS for densities 
n c < n < n s but softens it at higher densities. Their resulting maximum masses 
are in the range 2.2M < M < 3M Q . Our approach incorporates causality by 
reducing the sound speed smoothly towards the speed of light at high densities. 
Therefore, our maximum mass never exceeds that of the nuclear EoS of Akmal 
et al. [13]. In fact one may argue that at very high densities particles become 
relativistic and the sound speed should be even lower, v^ ~ c 2 /3. On the 
other hand, if it turns out that the QPOs are not from the innermost stable 
orbits and that even accreting neutron stars have small masses, say like the 
binary pulsars M^1.5M , this may indicate that heavier neutron stars are 
not stable. Therefore, the EoS is soft at high densities 5^,0.4 or that a phase 
transition occurs at a few times nuclear matter densities. For the nuclear to 
quark matter transition this would require 5 1 / 4 ^,80MeVfm~ 3 for 5 = 0.2. For 
such small Bag parameters there is an appreciable quark and nuclear matter 
mixed phase in the neutron star interior but even in these extreme cases a 
pure quark matter core is not obtained for stable neutron star configurations. 








Fig. 30. Total mass M for various values of 6. See text for further details. 
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Finally, we end this subsection by listing in Table 2 the maximum values 
for masses, radii and moments of inertia for several of the equations of state 
discussed in Figs. 30-38. From this table we see that the pn EoS with the lowest 
value of 6 gives also the stiffest EoS, and thereby largest mass and smallest 
central density. Similarly, the largest value for the bag constant results also 
in the stiffest EoS. In connection with the discussion of QPO's, it is worth 
pointing out that in Kerr space the relation between the Keplerian orbital 



frequency vk and the mass of the star is 2.198M (z/^kHz)" 



0.748J) 



with j = Iu/M a dimensionless measure of the angular momentum of the 
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Fig. 32. Total mass M for various values of the bag parameter B for the mixed 
phase EoS. For comparison we include also the results from the pn-matter EoS for 
/3-stable with 5 = 0.2. 
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Fig. 33. Mass-radius relation for various values of the bag parameter B for the mixed 
phase EoS. For comparison we include also the results from the pn-matter EoS for 
/3-stable with 5 = 0.2. 

staiQ. Following Ref. [139] and inserting the 1171 Hz QPO from 4U 1636-536, 
a rotational frequency u/2ir = 272 Hz and an assumed moment of inertia of 
~ lOOM km 2 results in a mass of 2.O2M and a radius of 9.6 ± 0.6 km. From 
the above table, we see that these results are fairly close to those which we 
get for the pure pn EoS with 6 = 0.2, i.e. for the 6 value which gave the best 



Recall that in all equations G = c = H = 1. 
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Fig. 34. Total mass M for various values of the bag parameter B for the Maxwell 
contructed EoS. For B 1 ' = 100 MeVfm" 3 , the purepn phase ends at 0.58 fm -3 and 
the pure quark phase starts at 0.67 fm" 3 . For B 1 ' 4 = 150 MeVfm" 3 , the numbers 
are 0.92 fm" 3 and 1.215 fm" 3 , while the corresponding numbers for B 1 ' 4 = 200 
MeVfm" 3 are 1.04 and 1.57 fm" 3 . In the density regions where the two phases 
coexist, the pressure is constant, a fact reflected in the constant value of the neutron 
star mass. 
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Fig. 35. Radius as function of central density n c for various values of the bag pa- 
rameter B for the Maxwell constructed EoS. For comparison we include also the 
results from the pn-matter EoS for /3-stable with 5 = 0.2. The region where R is 
constant reflects the density region where the pressure is constant in the Maxwell 
construction. See also the figure caption to the previous figure. 
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Fig. 36. Rotational mass M and gravitational mass for the pure pn EoS with 5 = 0.2 
and equations of state based on the mixed phase construction for different values 
of B. 
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Fig. 37. Moment of Inertia / in units of M0km as function of Mq for the pure pn 
EoS with 5 = 0.2 and for the mixed phase construction with B 1 ' 4 = 200 MeVfcm 3 . 

fit to the EoS of Akmal et al. [13]. Thus, if QPO's occur near the innermost 
stable orbits, then neutron star masses are M ~ 2.2M Q . This constrains the 
nuclear EoS including causality in a smooth way to only the stiffest ones - 
specifically 8^0.2 in the EoS of Eq. (49). Phase transitions in cores of neutron 
stars softens the EoS and strong transitions can therefore be ruled out except 
at very high densities n^,5n . On the other hand, if it turns out that the QPO 
are not from the innermost stable orbits and that even accreting neutron stars 
have small masses, say like the binary pulsars M^1.5M & , this indicates that 
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Fig. 38. Moment of Inertia I in units of M^km 2 as function of central density n c for 
the pn EoS with S = 0.2 and with quark degrees of freedom with B 1 ' = 200 
MeVfm -3 for the mixed phase and Maxwell constructions. Note well that the 
Maxwell construction yields a constant I since the pressure is constant in this case 
in the density region from 1.04 to 1.57 fm~ 3 . 



Table 2 

Maximum gravitational mass in Mq and the corresponding radius R, in units of 
km, for the given central density n c , in units of fm -3 , for various equations of state. 
The maximum moment of Inertia I, in units of M Q km 2 , is also listed. Note that 
this occurs for another central density than listed below. All results are for /3-stable 
matter and rotational corrections have not been included in the total mass. 
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Max mass 
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1.58 
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B 1 ^ = 200 MeVfm- 3 1.84 

B l / A = 150 MeVfm" 3 1.77 

B l / A = 100 MeVfm' 3 1.71 



91.5 11.6 1.1 
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heavier neutron stars are not stable. Therefore the EoS must be soft at high 
densities, i.e. 6^,0.4 or that a phase transition occurs at a few times nuclear 
matter densities. 

The role of phase transitions and its possible link with observation will be 
discussed in the following subsections. 



5.4 Maximum masses 



A stellar object of mass ~ 1.4M can either be an ordinary star of type F, a 
white dwarf, a neutron star, a black hole, or possibly a quark star, see Fig. 5.4. 
As shown by Fechner and Joss [9] a second branch of quark stars are possible 
for certain equation of states (EoS) - specifically for some parameter values of 
the Bag constant. The stars considered were hybrid stars consisting of a quark 
matter core with a mantle of nuclear matter around. A double maximum mass 
in the mass-density plot for neutron and mixed phase quark stars have also 
been found in [147] (see also [148]). 

The occurrence of a second maximum is a curious phenomenon that occurs 
under specific conditions that can be quantified. For that purpose we first con- 
sider a simple model that can be solved analytically, namely an EoS consisting 
of two incompressible fluids with a first order phase transition between energy 
density E\ and e 2 (£1 < £2) coexisting at a pressure Pq. We shall also first 
ignore effects of general relativity, i.e. take the Newtonian limit. 

The mass function, M(r) = An Jq e(r')r' 2 dr', is very simple in the Newtonian 
limit and the boundary condition M(R) = M relates the star mass M and 
radius R to the radius of the dense core, Rq, as 

M(R) = y (ei/2 3 + (e 2 - £i)#o) , (126) 

where Rq is the radius of the dense core. From Newton's equation for hydro- 
static equilibrium 



dP M(r 

dr r 



s^. (127) 



the pressure is easily obtained. From the boundary condition P(R) = we 
obtain 

F„ = | £ f(« 2 + (2|-3K-(|-l)|). (128) 

When the dense core is small, Rq -C R, the R^ terms in Eqs. (126,) and (128) 
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Fig. 39. Mass vs. central density for the Bethe-Johnson polytrope (T = 2.54) [149] 
with first order phase transitions at e\ = 3.2eo to e 2 - A region of instability occurs 
when Eq. (129) is fulfilled. 

can be ignored. From Eq. (128) we therefore observe that when 2e 2 > 3ei, 
the radius R of the star decreases with increasing size of the dense core Rq. 
Correspondingly, its mass M(R) of Eq. (126) decreases. In other words, as the 
average density of the star increases, its mass decreases and a stability analysis 
reveals that the star is unstable. It will contract until Rq is comparable to R 
such that the R^ terms in Eqs. (126) and (128) stabilizes the star and its mass 
again increases with increasing size. We can therefore conclude that when 



e 2 > - £l 



> instability region, 
a second maximum mass appears. 



(129) 



Another case, that can be solved analytically, is the T = 2 polytropic EoS, 
P oc e 2 . In this case the Newtonian version of the TOV equation is equivalent 
to the Schrodinger equation in a square well (or the Klein-Gordon equation) 
with solution e(r) oc sin(irr/R)/r. Including a first order phase transition 
leads to a phase shift in the corresponding sine solution for the outer mantle 
at densities e < S\. Curiously, one finds exactly the same instability criteria as 
Eq. (129). However, varying the polytropic index around T = 2 does change 
this condition slightly. Also, including general relativity affects the instability 
condition of Eq. (129) when e 2 is close to the maximum central density where 



•SO 



the star becomes unstable with respect to a collapse to a black hole, see Fig. 
5.4. A second order phase transition as, e.g. the mixed nuclear and quark 
matter phase, can also lead to instabilities and second maximum masses when 
the EoS is sufficiently softened [147]. 



5.5 Phase transitions in rotating neutron stars 



During the last years, and as discussed in the preceeding sections, interesting 
phase transitions in nuclear matter to quark matter, mixed phases of quark 
and nuclear matter [5,10], kaon [11] or pion condensates [13], neutron and 
proton superfluidity [8], hyperonic matter [14,16,17] crystalline nuclear mat- 
ter, magnetized matter, etc., have been considered. Recently, Glendenning et 
al. [6] have considered rapidly rotating neutron stars and what happens as 
they slow down when the decreasing centrifugal force leads to increasing core 
pressures. They find that a drastic softening of the equation of state, e.g. by a 
phase transition to quark matter, can lead to a sudden contraction of the neu- 
tron star at a critical angular velocity and shows up in a backbending moment 
of inertia as function of frequency. Here we consider another interesting phe- 
nomenon namely how the star and in particular its moment of inertia behaves 
near the critical angular velocity where the core pressure just exceeds that 
needed to make a phase transition. We calculate the moment of inertia, angu- 
lar velocities, braking index, etc. near the critical angular velocity and discuss 
observational consequences for first and second order phase transitions. 

Here we will make the standard approximation of slowly rotating stars, i.e. 
the rotational angular velocity is Q 2 <C M/R 3 . For neutron stars with mass 
M = 1.4M and radius R ~ 10 km their period should thus be larger than 
a few milliseconds, a fact which applies to all measured pulsars insofar. The 
general relativistic equations for slowly rotating stars were presented by Hartle 
[141] and reviewed in subsection 5.2. Hartle's equations are quite elaborate to 
solve as they consist of six coupled differential equations as compared to the 
single Tolman-Oppenheimer-Volkoff equation in Eq. (102) in the non-rotating 
case. In order to be able to analytically extract the qualitative behavior near 
the critical angular velocity Qq, where a phase transition occurs in the center, 
we will first solve the Newtonian equations for a simple equation of state. This 
will allow us to make general predictions on properties of rotating neutrons 
stars when phase transitions occur in the interior of a star. The corrections 
from general relativity are typically of order M/R ~10-20% for neutron stars 
of mass M ~ 1.4M . The extracted analytical properties of a rotating star 
are then checked below by actually solving Hartle's equations numerically for 
a realistic equation of state. 

The simple Newtonian equation of motion expresses the balance between the 



pressure gradient and the gravitational and centrifugal forces 

VP = -£(VV + fixfixr), (130) 

Here, V(r) is the gravitational potential for the deformed star and e the energy 
(~mass) density. We assume that friction in the (nonsuperfluid) matter insures 
that the star is uniformly rotating. Since cold neutron stars are barotropes, 
i.e. the pressure is a function of density, the pressure, density and effective 
gravitational potential, $ = V — |(J7 x r) 2 , are all constants on the same 
isobaric surfaces for a uniformly rotating star [141]. We denote these surfaces 
by the effective radius, a, and for slowly rotating stars it is related to the 
distance r from the center and the polar angle 9 from the rotation axis along 
Q by [141] 

r(a,0) = a[l-e(a)P 2 (cos0)], (131) 

where P2(cos8) is the 2nd Legendre polynomial and e(a) is the deformation 
of the star from spherical symmetry. 

Inserting Eq. (131) in Eq. (130) one obtains for small deformations [141] the 
I = Newtonian hydrostatic equation 

e da a z 3 

where M(a) = 4-7T J a e{a')a' 2 da' is the mass contained inside the mean radius 
a. The factor 2/3 in the centrifugal force arises because it only acts in two of 
the three directions. The equation (/ = 2) for the deformation e(o) is [150]. 

( l de \ *^~n ,- f ^rJe(a') _ 5 o2 



., , . M(a) - 4tt / p( a ')—^da' = -Q 2 . (133) 

a z da I J da' 3 



a 



The deformation generally increases with decreasing density, i.e. the star is 
more deformed in its outer layers. 

In order to discuss the qualitative behavior near critical angular velocities we 
first consider a simple EoS with phase transitions for which Eq. (132) can be 
solved analytically namely that of two incompressible fluids with a first order 
phase transition between energy density S\ and £2 (^i < £2) coexisting at a 
pressure Pq. The mass function M(a) is very simple in the Newtonian limit 
and the boundary condition M(R) = M relates the star radius R to the radius 
of the dense core, i? , as 

i?=(/2 3 -(|-l)^) V3 , (134) 



where R = (3M /Airei) 1 ^ 3 is the star radius in the absence of a dense core. 
Solving Eq. (132) gives the pressure 

P{a) = P + \{Rl - a 2 )e 2 (^e 2 - j^ 2 ), (135) 

for < a < R and 

1 4-ir 2 

P(a)=Po + -(R 2 -a*)e 1 (-e 1 --tf) 

-^( £2 _ £l ) £l (i_^), (136) 

3 a 

for R < a < R [7]. The boundary condition at the surface P(R) = in Eq. 
(136) gives 

2 1 <-> / 



lo z = = 1-2 



2ne 1 \Ane\R 2 

+ (| " 1)§(1 - ^o/i?)) (1 - Rl/R 2 Y l - (137) 

The phase transition occurs right at the center when Rq = corresponding to 
the critical angular velocity Qq = uj^Itiei where [7] 

^o = l-2^-. (138) 

Generally, for any EoS the critical angular velocity depends on Pq, M, and S\ 

but not on 62- 

For angular velocities just below uo$ very little of the high density phase exists 
and Rq -C R. Expanding (137) we obtain 



Rn 



2 _,.,2 

(139) 



UJq — UJ 



R \3- 2e 2 je x - u\ ' 

For uj > u>o the dense phase disappears and Rq = 0. Generaly, one can inter- 
prete Rq as an order parameter in analogy to, e.g. magnetization, the BCS 
gap, or the Higgs field in the standard model, however, as function of angu- 
lar velocity instead of temperature. Note, that for large density differences, 
£2/^1 > (3 — 0Oq)/2, Eq. (139) is not valid. This is related to an instability (see 
Eq. (129) and will be discussed in the following subsection. 

The corresponding moment of inertia is for R ^C R 

4tt / \ 2 

I=-(e 2 Rl + e 1 (R 5 -Rl))(l + -e) 




(140) 



where we used that the deformation from Eq. (133) is e = (5/4) u 2 in the 
low density phase [150]. However, for the qualitative behavior near Q only 
the contraction of the star radius R with the appearance of the dense core 
Ro is important whereas the deformations can be ignored. The contraction 
is responsible for the term in the moment of inertia and is proportional to 
Rq oc (to 2 — u) 3 / 2 near the critical angular velocity. Consequently, the derivative 
dl/du 2 displays the same non-analytic square root dependence as i? (see Eq. 
(139)). 

Latent heat is generated in the phase transition can be ignored because of 
rapid neutrino cooling which will be even faster than in supernova explosions. 
Thus temperatures will drop below ~ lMeV in seconds. Such temperatures 
are negligible compared to typical Fermi energies of nucleons or quarks and 
the timescales are also much smaller than to. 

Let us subsequently consider a more realistic EoS for dense nuclear matter at 
high densities such as the Bethe- Johnson EoS [149]. At high densities it can 
be approximated by a polytropic relation between the pressure and energy 
density: P = i^e: 2 ' 54 , where K\ = 0.021e^ L54 and Eq = m„0.15fm~ 3 is normal 
nuclear matter mass density. As we are only interested in the dense core we will 
for simplicity employ this Bethe- Johnson polytrope (BJP) EoS. The central 
density of a non-rotating 1.4M Q mass neutron star with the BJP EoS is ~ 
3.4e - Furthermore, we assume that a first order phase transition occurs at 
density E\ = 3.2e to a high density phase of density e<i = 4e with a similar 
polytropic EoS P = f^ 2 ' 54 - From the Maxwell construction the pressure is 
the same at the interface, Pq, which determines K.2 = K\{e\/ 'e^) 2 ' 54 '. We now 
generalize Eq. (132) by including effects of general relativity. From Einstein's 
field equations for the metric we obtain from the I = part 

1 iP _ M („) + WP + 2 



e + Pda a 2 (l-2M(a)/a) 3 

where m(a) = 4ti j a e{a')a' 2 da' . In the centrifugal force term we have ignored 
frame dragging and other corrections of order Q 2 M/R ~ 0.1f2 2 for simplicity 
and since they have only minor effects in our case. By expanding the pres- 
sure, mass function and gravitational potential in the difference between the 
rotating and non-rotating case, Eq. (141) reduces to the / = part of Hartle's 
equations (cf. Eq. (100) in [141].) Note also that Hartle's full equations cannot 
be used in our case because the first order phase transition causes disconti- 
nuities in densities so that changes are not small locally. This shows up, for 
example, in the divergent thermodynamic derivate de/dP. 



The rotating version of the Tolman-Oppenheimer-Volkoff equation (141) is 
now solved for a rotating neutron star of mass M = 1.4M with the BJP EoS 
including a first order phase transition. In Fig. 40 we show the central den- 
sity, moment of inertia, braking index, star radius and radius of the interface 
(Rq) as function of the scaled angular velocity. It is important to note that 




i?o oc A/ft 2 , — ^ 2 f° r angular velocities just below the critical value fto- The 
qualitative behavior of the neutron star with the BJP EoS and a first order 
phase transition is the same as for our simple analytic example of two incom- 
pressible fluids examined above. Generally, it is the finite density difference 
between the phases that is important and leads to a term in the moment of 
inertia proportional to (ft 2 , — ft 2 ) 3 / 2 as in Eq. (140). 

The moment of inertia increases with angular velocity. Generally, for a first 
order phase transition we find for ft,$,fto (see also Eq. (140) and Fig. (40)) 



(142) 



For the two incompressible fluids with momentum of inertia given by Eq. 
(140), the small expansion parameters are c\ = uJq and c 2 = (5/2)cJq(£ 2 /£i — 
l)/(3 — 2e 2 /e 1 — Wq) 3//2 ; for ft > ft the c 2 term is absent. For the BJP we find 
from Fig. (1) that c 2 — 0.07 ~ 2. 2a; 3 ,. Generally, we find that the coefficient 
c 2 is proportional to the density difference between the two coexisting phases 
and to the critical angular velocity to the third power, c 2 ~ {&ijs\ — 1 Vq. The 
scaled critical angular velocity ujq can at most reach unity for submillisecond 
pulsars. 

To make contact with observation we consider the temporal behavior of an- 
gular velocities of pulsars. The pulsars slow down at a rate given by the loss 
of rotational energy which we shall assume is proportional to the rotational 
angular velocity to some power (for dipole radiation n — 3) 

With the moment of inertia given by Eq. (142) the angular velocity will then 
decrease with time as 



n en 



n — 1 



ft I 








n-l)t \ V t 




' '• ' (^)2/("-l) + .... I , (144) 



for t > t . Here, the time after formation of the pulsar is, using Eq. (143), 
related to the angular velocity as t ~ t (fto/ft) n_1 and t = lQ/({n — l)CVtQ~ l ) 
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for n > 1, is the critical time where a phase transition occurs in the center. 
For earlier times t < t there is no dense core and Eq. (144) applies when 
setting c 2 = The critical angular velocity is Q = oj^Iixex — 6kHz for 
the BJP EoS, i.e. comparable to a millisecond binary pulsar. Applying these 
numbers to, for example, the Crab pulsar we find that it would have been 
spinning with critical angular velocity approximately a decade after the Crab 
supernova explosion, i.e. to ~ 10 years for the Crab. Generally, to oc Qo~ n and 
the timescale for the transients in Ci as given by Eq. (13) may be months or 
centuries. In any case it would not require continuous monitoring which would 
help a dedicated observational program. 

The braking index depends on the second derivative I" = dI/d 2 Q of the 
moment of inertia and thus diverges (see Fig. (1)) as Q approaches Qq from 
below 

n ( fi ) = ™ „ n - 2 Cl % + c 2 — r £22L= . (145) 

For Q > Qq the term with c 2 is absent. The observational braking index n(Q) 
should be distinguished from the theoretical exponent n appearing in Eq. (12). 
Although the results in Eqs. (13) and (14) were derived for the pulsar slow 
down assumed in Eq. (12) both Q and n(Q) will generally display the y/t — t 
behavior for t > to as long as the rotational energy loss is a smooth function 
of Q. The singular behavior will, however, be smeared on the pulsar glitch 
"healing" time which in the case of the Crab pulsar is of order weeks only. 

We now discuss possible phase transitions in interiors of neutron stars. The 
quark and nuclear matter mixed phase described in [5] has continuous pres- 
sures and densities. There are no first order phase transitions but at most two 
second order phase transitions. Namely, at a lower density, where quark matter 
first appears in nuclear matter, and at a very high density (if gravitationally 
stable), where all nucleons are finally dissolved into quark matter. In second- 
order phase transitions the pressure is a continuous function of density and we 
find a continuous braking index. This mixed phase does, however, not include 
local surface and Coulomb energies of the quark and nuclear matter struc- 
tures. As shown in [10,135] there can be an appreciable surface and Coulomb 
energy associated with forming these structures and if the interface tension 
between quark and nuclear matter is too large, the mixed phase is not favored 
energetically. The neutron star will then have a core of pure quark matter with 
a mantle of nuclear matter surrounding it and the two phases are coexisting 
by a first order phase transition. For a small or moderate interface tension the 
quarks are confined in droplet, rod- and plate-like structures [10,135] as found 
in the inner crust of neutron stars [33]. Due to the finite Coulomb and surface 
energies associated with forming these structures, the transitions change from 
second order to first order at each topological change in structure. If a Kaon 
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condensate appears it may also have such structures [91]. Pion condensates 
[153], crystalline nuclear matter [12,13], hyperonic or magnetized matter, etc. 
may provide other first order phase transitions. 

If a neutron star cools continuously, the temperature will decrease with time 
and the phase transition boundary will move inwards. The two phases could, 
e.g. be quark-gluon/nuclear matter or a melted/solid phase. In the latter case 
the size of the hot (melted) matter in the core is slowly reduced as the tem- 
perature drops freezing the fluid into the solid mantle. Melting temperatures 
have been estimated in [33,133] for the crust and in [10] for the quark matter 
mixed phase. When the very core freezes we have a similar situation as when 
the star slows down to the critical angular velocity, i.e. a first order phase 
transition occurs right at the center. Consequently, a similar behavior of mo- 
ment of inertia, angular velocities, braking index may occur as in Eqs. (142), 
(144) and (145) replacing Q(t) with T(t). 

Thus, if a first order phase transitions is present at central densities of neutron 
stars, it will show up in moments of inertia and consequently also in angular 
velocities in a characteristic way. For example, the slow down of the angular 



velocity has a characteristic behavior fi ~ ci\j 1 — t/t and the braking index 



diverges as n(f2) ~ c 2 /Jl — Vt 2 /VLq (see Eqs. (144) and (145)). The magnitude 
of the signal generally depends on the density difference between the two 
phases and the critical angular velocity ujq = Qo/^/2ne^ such that c 2 ~ (£2/£i — 
1)u;q. The observational consequences depend very much on the critical angular 
velocity f2 , which depends on the equation of state employed, at which density 
the phase transition occurs and the mass of the neutron star. 

We encourage a dedicated search for the characteristic transients discussed 
above. As the pulsar slows down over a million years, its central densities 
spans a wide range of order Iuq (see Fig. 40). As we are interested in time 
scales of years, we must instead study all ~ 1000 pulsars available. By study- 
ing the corresponding range of angular velocities for the sample of different 
star masses, the chance for encountering a critical angular velocity increases. 
Eventually, one may be able to cover the full range of central densities and find 
all first order phase transitions up to a certain size determined by the experi- 
mental resolution. Since the size of the signal scales with Qq the transition may 
be best observed in rapidly rotating pulsars such as binary pulsars or pulsars 
recently formed in supernova explosion and which are rapidly slowing down. 
Carefully monitoring such pulsars may reveal the characteristic behavior of 
the angular velocity or braking index as described above which is a signal of 
a first order phase transition in dense matter. 
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Fig. 40. Central density (in units of £o), radii of the neutron star R and its dense 
core Rq, moment of inertia, its derivative I' /I = dl/dco 2 /I and the braking index 
are shown as function of the scaled angular velocity io 2 = Q 2 /(2irei). The rotating 
neutron star has mass IAMq and a Bethe- Johnson like polytropic equation of state 
with a first order phase transition taking place at density e± = 3.2eo to £2 = 4eo- 

5. 6 Core quakes and glitches 



The glitches observed in the Crab, Vela, and a few other pulsars are probably 
due to quakes occurring in solid structures such as the crust, superfluid vortices 
or possibly the quark matter lattice in the core [135] . As the rotating neutron 
star gradually slows down and becomes less deformed, the rigid component 
is strained and eventually cracks/quakes and changes its structure towards 
being more spherical. 



The moment of inertia of the rigid component, I c , decreases abruptly and its 
rotation and pulsar frequency increases due to angular momentum conserva- 
tion resulting in a glitch. The observed glitches are very small AQ/Q ~ 10 -8 . 
The two components slowly relaxate to a common rotational frequency on a 
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timescale of days (healing time) due to superfluidity of the other component 
(the neutron liquid). The healing parameter Q = I c /Itot measured in glitches 
reveals that for the Vela and Crab pulsar about ~3% and ~96% of the moment 
of inertia is in the rigid component respectively. 

If the crust were the only rigid component the Vela neutron star should be 
almost all crust. This would require that the Vela is a very light neutron star - 
much smaller than the observed ones which all are compatible with ~ 1.4M . 
If we by the lattice component include not only the solid crust but also the 
protons in nuclear matter (NM) (which is locked to the crust due to magnetic 
fields), superfluid vortices pinned to the crust [154] and the solid QM mixed 
phase 

I C = I crust + Ip + Isv + IQM , (146) 

we can better explain the large I c for the Crab. The moment of inertia of 
the mixed phase is sensitive to the EoS's used. For example, for a quadratic 
NM EoS [10] decreasing the Bag constant from 110 to 95 MeVfm -3 increases 
I d I total from ~ 20% to ~ 70% for a 1.4M Q neutron star - not including possible 
vortex pinning. The structures in the mixed phase would exhibit anisotropic 
elastic properties, being rigid to some shear strains but not others in much 
the same way as liquid crystals. Therefore the whole mixed phase might not 
be rigid. 

The energy released in glitches every few years are too large to be stored in the 
crust only. The recurrence time for large quakes, t c , is inversely proportional 
to the strain energy [154], which again is proportional to the lattice density 
and the Coulomb energy 

^oc^— ■ (147) 

a 6 a 

Since the lattice distance a is smaller for the quark matter droplets and their 
charge larger than for atoms in the crust, the recurrence time is shorter in 
better agreement with measurements of large glitches. 

Detecting core and crust quakes separately or other signs of three components 
in glitches, indicating the existence of a crust, superfluid neutrons and a solid 
core, would support the idea of the mixed quark and nuclear matter mixed 
phase. However, magnetic field attenuation is expected to be small in neutron 
stars and therefore magnetic fields penetrate through the core. Thus the crust 
and core lattices as well as the proton liquid should be strongly coupled and 
glitch simultaneously. 
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5. 7 Backbending and giant glitches 



In [6] the moment of inertia is found to "backbend" as function of angular 
velocity. The moment of inertia of some deformed nuclei [151,152] may also 
backbend when the coriolis force exceeds the pairing force breaking the pair- 
ing whereby the nucleus reverts from partial superfluidity to a rigid rotor. 
However, in the limit of large nuclear mass number such backbending would 
disappear. Instead pairing may lead to superfluidity in bulk [154]. A backbend- 
ing phenomenon in neutron stars, that appears to be similar to backbending 
in nuclei, can occur in neutron stars although the physics behind is entirely 
different. If we soften the EoS significantly at a density near the central den- 
sity of the neutron star, a non-rotating neutron star can have most of its core 
at high densities where the soft EoS determines the profile. A rapidly rotating 
star may instead have lower central densities only probing the hard part of the 
EoS. Thus the star may at a certain angular velocity revert from the dense 
phase to a more dilute one and at the same time change its structure and 
moment of inertia discontinuosly. Such a drastic change in moment of inertia 
at some angular velocity will cause a giant glitch as found in [6]. 

The phenomenon of neutron star backbending is related to the double max- 
imum mass for a neutron and quark star respectively as shown in Fig. (5.4). 
The instabilities given in Eq. (129) are also evident when rotation is included 
as seen from Eq. (139) 

£2 > — ^i(3 — uOq) =>- Giant glitch when 00 ~ lu . (148) 

The Bethe- Johnson EoS discussed in subsection 5.4 also has this discontinuity 
in the moment of inertia when p 2 /Pi~3/2. The neutron star may continue to 
slow down in its unstable structure, i.e. "super-rotate" , before reverting to its 
stable configuration with a dense core. As for the instabilities of Eq. (129) this 
condition changes slightly for a more general EoS and when general relativity 
is included. Neutron stars with a mixed phase do not have a first order phase 
transition but may soften their EoS sufficiently that a similar phenomenon 
occurs [6]. We emphasize, however, that the discontinuous jump in moment of 
inertia is due to the drastic and sudden softening of the EoS near the central 
density of neutron stars. It is not important whether it is a phase transition 
or another phenomenon that causes the softening. 



5.8 Cooling and temperature measurements 



The thermal evolution of a neutron star may provide information about the 
interiors of the star, and in recent years much effort has been devoted to mea- 
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Table 3 

Luminosities, L, and spin-down ages, r, of pulsars 



Pulsar 


Name 


logr [yr] 


logL [erg/s] 


1706-44 




4.25 


32.8 ± 0.7 [156] 


1823-13 




4.50 


33.2 ± 0.6 [157] 


2334+61 




4.61 


33.1 + 0.4 [158] 


0531+21 


Crab 


3.09 


33.9 ± 0.2 [159] 


1509-58 


SNR MSH 15-52 


3.19 


33.6 ± 0.4 [160] 


0540-69 




3.22 


36.2 + 0.2 [161] 


1951+32 


SNR CTB 80 


5.02 


33.8 ± 0.5 [162] 


1929+10 




6.49 


28.9 + 0.5 [163,155] 


0950+08 




7.24 


29.6 ± 1.0 [164] 


J0437-47 




8.88 


30.6 ± 0.4 [165] 


0833-45 


Vela 


4.05 


32.9 ± 0.2 [166] 


0656+14 




5.04 


32.6 ± 0.3 [167] 


0630+18 


Geminga 


5.51 


31.8 + 0.4 [168] 


1055-52 




5.73 


33.0 ± 0.6 [169] 



suring neutron star temperatures, especially with the Einstein Observatory 
and ROSAT, see e.g. Ref. [155]. Neutron stars are born with interior tem- 
peratures of the order 20-50 MeV, but cool rapidly via neutrino emission to 
temperatures of the order 1 MeV within minutes. The only information on 
neutron star temperatures stems from surface temperatures, which typically 
are of the order 10 6 K for about 10 5 yr, observed in X-ray or UV bands. 
However, the thermal radiation from a neutron star has yet to be identified 
unambiguously. Most observations are for pulsars, and it is unclear how much 
of the observed radiation is due to pulsar phenomena, to a synchotron-emitting 
nebula or to the neutron star itself. Surface temperatures of neutron stars have 
been measured in a few cases or upper limits have been set. Table 3 collects 
some of these. 

The cooling history of the star, and energy loss mechanisms from the interior 
are thus to be determined through various theoretical models. The generally 
accepted picture is that the long-term cooling of a neutron star consists of two 
periods: a neutrino cooling epoch which can last until 10 6 yr and a photon 
cooling period. If we now assume that the main cooling mechanism in the 
early life of a neutron star is believed to go through neutrino emissions in the 
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core, the most powerful energy losses are expected to be given by the so-called 
direct Urea mechanism 

n — > p + l + ui, p + l—tn + vi, (149) 

as rediscussed recently by several authors [170-172]. The label / refers to the 
leptons considered here, electrons and muons. However, in order to fullfil the 
momentum conservation k F < k F + kp and energy conservation requirements, 
the process can only start at densities n several times nuclear matter saturation 
density no = 0.16 fm~ 3 , see e.g. Fig. 18, where the proton fraction exceeds 
x p >0.14, see e.g. [46,170-172]. 

Thus, for long time the dominant processes for neutrino emission have been 
the so-called modified Urea processes first discussed by Chiu and Salpeter 
[173], in which the two reactions 

n + n — ► p + n + I +Vi, p + n + l^n + n + vi, (150) 

occur in equal numbers. These reactions are just the usual processes of neutron 
/9-decay and electron and muon capture on protons of Eq. (149), with the 
addition of an extra bystander neutron. They produce neutrino-antineutrino 
pairs, but leave the composition of matter constant on average. Eq. (150) 
is referred to as the neutron branch of the modified Urea process. Another 
branch is the proton branch 

n + p — > p + p + l +Vi, p + p + l-^n + p + ui, (151) 

pointed out by Itoh and Tsuneto [174] and recently reanalyzed by Yakovlev 
and Levenfish [175]. The latter authors showed that this process is as effi- 
cient as Eq. (150). In addition one also has the possibility of neutrino-pair 
bremsstrahlung. These processes form the basis for what is normally called 
the standard cooling scenario, i.e. no direct Urea processes are allowed. 

A fast cooling scenario would involve the direct Urea process, similar direct 
processes with baryons more massive than the nucleon participating, such as 
isobars or hyperons [171,172], or neutrino emission from more exotic states 
like pion and kaon condensates [87,35,92] or quark matter [29,126]. Actually, 
Prakash et al. [171] showed that the hyperon direct Urea processes gave a 
considerable contribution to the emissivity, without invoking exotic states or 
the large proton fractions needed in Eq. (149). 

If we now consider the neutron stars discussed in subsection 5.3, we note that 
a typical star with mass ~ 1.4 — 2.OM has central densities ranging from 
~ 0.5 — 1.2 fm -3 . Depending on the value of the bag constant, the mixed 
phase could start already at ~ 0.2 — 0.3 fm~ 3 . No pure quark phase was found 
with the bag-model for stable neutron star configurations. If one also recalls 
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that our /3-stable EoS allows for the direct Urea process at densities starting 
from 0.8 fm -3 , see Fig. 18, one clearly sees that, depending on the EoS and the 
adopted model for dense matter, there is a considerable model dependence. 
To give an example, the cooling could be strongly influenced by the mixed 
phase. This could come about because nuclear matter in the droplet phase 
has a higher proton concentration than bulk neutral nuclear matter and this 
could make it easier to attain the threshold condition for the nucleon direct 
Urea process. Another possibility is that the presence of the spatial structure 
of the droplet phase might allow processes to occur which would be forbidden 
in a translationally invariant system. Also the mere presence of quark matter 
can lead to fast cooling [126] when a s ^ 0. All these mechanisms would lead 
to faster cooling. 

However, in order to compare with observation, the structure of the star has 
to be computed in detail. In particular the possible presence of superfluidity 
in the interior has to be considered. The superfluid would suppress the v, v 
emissivity and would allow for reheating through friction with the crust. In the 
analysis of Page as well [176], it is hard to discriminate between fast and slow 
cooling scenarios, though in both cases agreement with the observed temper- 
ature of Geminga is obtained if baryon pairing is present in most, if not all of 
the core of the star. The recent analyses of Schaab et al. [177,178] also seem to 
confirm the importance of superfluidity in the interior of stars. There are also 
indications [179,180] that temperatures of young (~ 10 4 years old) neutron 
stars lie below that obtained through the so-called modified Urea processes. 
One has also to note [176] that the modified Urea processes are weakly depen- 
dent on the mass of the star, i.e. on the central density, while faster cooling 
mechanisms like the above direct Urea processes are in general strongly depen- 
dent on it. Thus, the detection of two coeval stars, whose temperatures differ 
by a factor of the order of 2 or larger would allow to distinguish between tra- 
ditional cooling scenarios, like those discussed by Page [176], and more exotic 
ones. Such a huge variation in the temperature of coeval stars could indicate 
the presence of a threshold in the cooling mechanism, triggered by the density 
of the star. 

In Fig. 5.8 we display results from various cooling calculations by Schaab et 
al. [178]. All calculations employ a superfluidity scenario which is described in 
section 2.7. This corresponds to the label 4 in the above figure, while the letters 
A,B,C,D represent different equations of state employed to calculate the mass 
of the star. Label A corresponds to the non-relativistic equations of state from 
Wiringa et al. [84], and is similar to that of Akmal et al. [13] in section 2, and 
non-relativistic equations of state based on two-body interactions only. The 
models encompassed by class A allow for the standard cooling scenario only 
at higher densities. These densities however, are beyond the central density of 
a 1.4M Q neutron star. Models B and C allow for faster cooling scenarios and 
the equations of state are based on relativistic mean field models. Typically, 




Fig. 41. Thermal evolution of a 1.4M neutron star with superfluidity. See text for 
further details. Taken from Ref. [178]. 

see also the discussion of relativistic effects and the proton fractions in Fig. 
16, these models allow for the direct Urea for nucleons at lower densities than 
the non-relativistic models. Pion and kaon condensantion could also lead to 
faster cooling scenarios. Model D allows for direct hyperon cooling as well but 
does not include hyperon pairing, which reduces strongly the direct Urea for 
hyperons, see Ref. [110]. 

The problem with most cooling calculations is that there is no consistent 
calculation of properties entering the neutrino emissivities within the frame- 
work of say one given EoS and many-body approach. Typically, see e.g. Refs. 
[177,178] and Fig. 5.8, the EoS is taken from one source, while the pairing gap 
is taken from another calculation, with even entirely different NN interactions 
or many-body approaches. In addition, the expressions for the emissivities of 
e.g. the modified Urea processes calculated by Friman and Maxwell [181] treat 
in a rather cavalier way the role of many-body correlations. Considering also 
the fact that other severe approximations are made, these expressions, which 
enter typically various cooling codes, could introduce errors at the level of 
orders of magnitude. 

Thus, our message is that, before one attempts at a cooling calculation, little 
can be learned unless the various neutrino emissivities are reevaluated within 
the framework of a given many-body scheme for dense matter. 
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We conclude this subsection with a demonstration of the role of superfluidity 
for the processes of Eqs. (150)-(151) at densities corresponding to the outer 
core of massive neutron stars or the core of not too massive neutron stars 
when we have a superfluid phase. Here we limit ourselves to study the role 
of superconducting protons in the core of the star employing the gap for pro- 
tons in the 1 5'o and effective masses from lowest-order Brueckner-Hartree-Fock 
calculations discussed in section 2 and subsections 2.2 and 2.7. The proton su- 
perconductivity reduces the energy losses considerably in the above reactions 
[175], and may have important consequences for the cooling of young neutron 
stars. The expressions for the processes in Eqs. (150)-(151) were derived by 
Friman and Maxwell [181] and read [175] 

Q n » 8.5 x 10 21 (^V (^i) (^-) 1/3 T*a n (3 n , (152) 

\m n J \m p J \n / 

in units of ergs cm _3 s _1 where Tg is the temperature in units of 10 9 K, and 
according to Friman and Maxwell a n describes the momentum transfer de- 
pendence of the squared matrix element in the Born approximation for the 
production rate in the neutron branch. Similarly, f3 n includes the non-Born cor- 
rections and corrections due to the nucleon-nucleon interaction not described 
by one-pion exchange. Friman and Maxwell [181] used a n ~ 1.13 at nuclear 
matter saturation density and (3 n = 0.68. In the results presented below, we 
will not include a n and f3 n . For the reaction of Eq. (150) with muons, one has 
to replace n e with n M and add a factor (1 + X) with X = k F /kp [182,183]. 
For the proton branch of Eq. (151) with electrons we have the approximate 
equation [175] 

Q p « 8.5 x 10 21 (<] 3 (!<) (!±) V3 T$a p (3 p F, (153) 

\m p ) \m n J \n J 

with F — I 1 — jpr ) 0, where = 1 if k F < 3k F + kp and zero elsewhere. 

Yakovlev and Levenfish [175] put a n = a p and /3 p = /3 n . We will, due to 
the uncertainty in the determination of these coefficients, omit them in our 
calculations of the reaction rates. With muons, the same changes as in Eq. 
(152) are made. 

The reaction rates for the Urea processes are reduced due to the supercon- 
ducting protons. Here we adopt the results from Yakovlev and Levenfish [175], 
their Eqs. (31) and (32) for the neutron branch of Eq. (152) and Eqs. (35) 
and (37) for the proton branch of Eq. (153). We single out proton singlet- 
superconductivity only, employing the approximation ksT c = A(0)/1.76, 
where Tq is the critical temperature and A(0) is the pairing gap at zero tem- 
perature discussed in section 2. The critical temperature is then used to obtain 
the temperature dependence of the corrections to the neutrino reaction rates 
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Fig. 42. Temperature dependence of neutrino energy loss rates in a neutron star 
core at a total baryonic density of 0.3 fm -3 with muons and electrons. Solid line 
represents Eq. (2) with electrons, dotted line is Eq. (2) with muons, dashed line is 
Eq. (3) with electrons while the dash-dotted line is the processes of Eq. (3) with 
muons. The corresponding results with no pairing are also shown. 

due to superconducting protons. To achieve that we employ Eq. (23) of [175]. 
We present results for the neutrino energy rates at a density n = 0.3 fm~ 3 . 
The critical temperature is Tq = 3.993 x 10 9 K, whereas without muons we 
have Tq = 4.375 x 10 9 . The implications for the final neutrino rates are shown 
in Fig. 42, where we show the results for the full case with both muons and 
electrons for the processes of Eqs. (150) and (151). In addition, we also display 
the results when there is no reduction due to superconducting protons. At the 
density considered, n = 0.3 fm -3 , we see that the processes of Eqs. (150) and 
(151) with muons are comparable in size to those of Eqs. (150) and (151) with 
electrons However, the proton pairing gap is still sizable at densities up to 
0.4 fm" 3 , and yields a significant suppression of the modified Urea processes 
discussed here, as seen in Fig. 42. We have omitted any discussion on neutron 
pairing in the 3 P2 state. For this channel we find the pairing gap to be rather 
small, see again the discussion in subsection 2.7, less than 0.1 MeV and close 
to that obtained in Ref. [102]. We expect therefore that the major reductions 
of the neutrino rates in the core come from superconducting protons in the 1 S'o 
state. The contribution from neutrino-pair bremsstrahlung in nucleon-nucleon 
collisions in the core is also, for most temperature ranges relevant for neutron 
stars, smaller than the contribution from modified Urea processes [175]. Pos- 
sible candidates are then direct Urea processes due to hyperons and isobars, 
as suggested in Ref. [110,171], or neutrino production through exotic states of 
matter, like kaon or pion condensation [35,87,92] or quark matter [29,126]. 

In conclusion, superfluidity reduces the neutrino emissivities considerably for 
the standard cooling model. However, before a firm conclusion from cooling 
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calculations can be obtained, one needs much more reliable estimates of various 
emissivity processes in dense matter. In addition, such processes should be 
calculated consistently within the same many-body approach and EoS model 
for the interior of a star. 



5. 9 Supernovae 



Neutron stars are born in type II or lb supernovae at a rate of 1-3 per century 
in our galaxy. Once the iron core in massive stars exceed the correspond- 
ing Chandrasekhar mass, it collapses adiabatically until it is stopped by the 
incompressibility of nuclear matter, bounces, creates shock waves, stalls by 
infalling matter, and presumably explodes by neutrinos blowing the infalling 
matter off after a few tens of seconds. 

Measurements of isotope abundances of various elements (r- and s-processes) 
give some insight in densities and temperatures during certain stages of the 
explosions. Neutrinos were also detected from SN1987A in the Large Magellan 
Cloud. 

Neutron stars are normally treated as an input to supernova calculations of 
the complicated transport processes of baryons, photons, leptons and neutri- 
noes going on during the explosion. Thus supernovae do not provide much 
information on details on neutron star structure presently. Recent supernova 
models include convection and spherical asymmetries (see, e.g. [184]). 

The softer the equation of state the denser the matter is compressed before 
it bounces and the deeper into the gravitational well the star has fallen. Also 
a softer EoS creates a more coherent shock wave that excites the matter less. 
The additional gravitational energy available can be transferred to neutrino 
generation which is believed to power the supernova explosion. Besides soft- 
ening the EoS a first order phase transition would store latent heat which also 
could affect supernovae. 

With luck we may observe a supernova nearby in the near future which pro- 
duces a rapidly rotating pulsar. Light curves and neutrino counts will test 
supernova and neutron star models. The rapid spin down may be exploited to 
test the structure and possible phase transitions in the cores of neutron stars 

[5-7]. 
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5.10 Gamma ray bursters 



In the late 1960's the Vela satellite was launched carrying a gamma ray de- 
tector in order to check the nuclear test ban treaty in space. It was a huge 
success because numerous gamma ray bursts (GRB) were observed [185] and 
(more importantly) they did not come from Earth. Later (when declassified) 
the russians confirmed the GRB. Numerous observations by the Pioneer Venus 
Orbiter, Compton Gamma Ray Observatory, Burst And Transient Source Ex- 
periment, and later gamma ray detectors now find bursts every day. The GRB 
do not repeat (except for a few soft gamma ray repeaters). Their duration 
varies from milliseconds to minutes. The short bursts imply a small source 
size £c ■ 1 ms~ 100km, which points towards neutron stars or black holes. 
The gamma rays have energy in the range 30 keV-2 MeV. The high energy 
gamma rays are above the threshold for 7 — »■ e + e~ photo production which 
implies that the radiation probably is relativisticaly expanding. 

The thousands of GRB observed show a high isotropy. Thus GRB cannot be 
produced in the very anisotropic discs of our solar system or galaxy or even an 
extended galactic halo. The final kill to such (and many other) models came 
from the Bepposax observations on may 8th, 1997, where a burst could be pin- 
pointed on the sky within an arcminut, which subsequently allowed ground 
based observations of an afterglow in optical and radio wavelengths. Fe-II and 
Mg-II absorption lines were found at high redshifts of z — 0.835. Assuming 
isotropic flux gives an energy output of -E 7 ^,10 51 ergs within seconds - much 
more powerful than supernovae. In comparison supernova produce ~ 10 51 erg's 
mainly in neutrinos whereas a long distance missile has (only) ~ 10 31 ergs of 
explosive energy. At the time of writing, two dozens of GRB with afterglow 
in X-ray, one dozen in optical and half a dozen in radio waves have been 
discovered. 

Merging neutron stars may be responsible for Gamma Ray Bursts [184]. Binary 
pulsars are rapidly spiralling inwards and will eventually merge and create a 
gigantic explosion and perhaps collapse to a black hole. From the number of 
binary pulsars and spiral rates, one estimates about one merger per million 
year per galaxy. With ~ 10 9 galaxies at cosmo logical distances z£l this gives 
about the observed rate of GRB. However, detailed calculations have problems 
with baryon contamination, i.e. the baryons ejected absorb photons in the 
relativistically expanding photosphere. Accreting black holes have also been 
suggested to act as beaming GRB. 

Sofar, the physics producing these GRB is not understood. The time scales 
and the enormous power output points towards neutron star or black hole 
objects. We hope to learn more about these objects and maybe someday they 
can be used to gain information on neutron stars. 
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6 Conclusions 



The aim of this work has been to give a survey of recent progresses in the 
construction of the equation of state for dense neutron star matter and possible 
implications for phase transitions inside neutron stars and the connections 
to neutron star observables. We will here try to recapitulate several of the 
arguments presented. 



6.1 Many-body approaches to the equation of state 



In section 2 we attempted at a review of the present status of microscopic 
many-body approaches to dense neutron star matter. Several features emerged: 

• Within non-relativistic lowest-order Brueckner theory (LOB), all the new 
phase-shift equivalent nucleon-nucleon potentials yield essentially similar 
equations of state up to densities of 3 — 4n for both pure neutron matter 
and /3-stable matter. Other properties like the symmetry energy and proton 
fractions do also show a similar quantitative agreement. 

• The inclusion of more complicated many-body terms at the two-body level 
does not alter this picture and even the recent summation of three-hole line 
diagrams of Baldo and co-workers [53,54] results in an EoS which is too close 
to LOB when a continuous choice is used for the single-particle energies in 
matter. 

• At densities up to nuclear matter saturation density n , basically all many- 
body approaches discussed here give very similar equations of state. 

• In symmetric nuclear matter the situation is however different, since there 
the nuclear tensor force is much more dominant due to presence of the 3 S\ 
and 3 Di partial waves. This has also consequences for the three-body in- 
teractions, arising from both effective and real three-body force terms. This 
is expected since the strong isospin T = channel is present in the three- 
body nnp and npp clusters. In neutron matter however, such clusters are in 
general small at densities up to n^. Similar behaviors are also seen in recent 
works on the energy of pure neutron drops [186] and larges-scale shell-model 
calculations of Sn isotopes including effective three-body interactions [187]. 
However, for PNM and densities greater than no differences do however oc- 
cur at densities when one introduces real three-body forces and/or includes 
relativistic corrections. In a similar way, relativistic BHF calculations, yield 
also significant corrections above saturation density. 

• Based on the microscopic calculation of Akmal et al. [13] , a simple parametriza- 
tion of the EoS was given in the previous subsection, where the energy per 
particle could be written as 

S = £ u h Sou'x 2 , 

1 + du 
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with So = 15.8MeV, 8 ~ 0.2, symmetry energy So = 32MeV and 7 = 0.6. 
Causality is insured by 5 > 0.13. 

In the case of non-nucleonic degrees of freedom we discussed also recent pro- 
gresses in the construction of hyperon-hyperon interactions [18,19] and ap- 
plications to neutron star matter [14]. Results for /3-stable matter with these 
recent hyperon-hyperon interactions are in progress [86]. In this work how- 
ever, we chose to focus on the Bag model in order to deal with non-nucleonic 
degrees of freedom. 

For studies of superfluidity in neutron star matter, one needs still a careful 
analysis of polarization effects, especially for the 1 S'o proton contaminant. For 
superfluid neutrons in the 3 P2 partial wave, one needs nucleon-nucleon in- 
teractions which are fitted up to at least 1 GeV in lab energy before firm 
conclusions can be reached. These points leave clearly a large uncertainty for 
the role of superfluidity in cooling studies. 



6.2 Phase transitions in neutron stars 



We have discussed a number of possible phase transitions in dense nuclear 
matter such as pion, kaon and hyperon kondensation, superfluidity and quark 
matter. We have specifically treated the nuclear to quark matter phase transi- 
tion and the possible mixed phase that can occur in more than one component 
systems and replace the standard Maxwell construction. The structure of the 
mixed phase is similar to the mixed phase in the inner crust of neutron stars 
of nuclei and the neutron gas. However, as the mixed quark and nuclear mat- 
ter phase can occur already at a few normal nuclear matter densities, it can 
soften the EoS in cores of neutron stars significantly and lower the maximum 
mass. A number of numerical calculations of rotating neutron stars with and 
without phase transitions were given. 

The calculated maximum masses were discussed with the observed ones and 
leave two natural options: 

• Case I: The large masses of the neutron stars in QPO 4U 1820-30 (M = 
2.3M & ), PSR J1012+5307 (M = 2.1 ± 0AM & ), Vela X-l (M = 1.9 ± 
O.IMq), and Cygnus X-2 (M = 1.8 ± 0.2Mq), turn out to be correct and 
are complemented by other neutron stars with masses around ~ 2M . 

As a consequence, the EoS of dense nuclear matter is severely restricted 
and only the stiffest EoS allowed by causality are allowed (i.e. 5 ~ 0.2). Also 
any significant phase transition can be excluded at densities below ^5no- 

That the radio binary pulsars all have masses around 1.4M is then 
probably due to the formation mechanism in supernovae and related to 
the Chandrasekhar mass M c ^ ~ 1.4M of white dwarfs. Neutron stars can 
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subsequently acquire larger masses by accretion. 
• Case II: The heavy neutron stars proves erroneous by more detailed obser- 
vations and only the canonical ~ 1.4M Q masses are found. 

If even in accreting neutron stars does not produce neutron stars heavier 
than say ^1.5M , this indicates that heavier neutron stars simply are not 
stable whic in turn implies a soft EoS, either 5^,0.4 or a significant phase 
transition. 

Either way the result is important information on neutrons stars, the EoS of 
dense hadronic matter and possible phase transitions. We impatiently await 
future observations and determinations of neutron stars, that will answer these 
questions. 
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